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Abstract 


In this thesis, we study the d —orthogonal polynomials and the d— orthogonality in the 
Sobolev sens. We present an algebraic theory of classical d—orthogonal polynomials and 
we want to fill in some gaps. We broaden and close some inclusions that exist and are 
known perhaps as consequences. Several characterizations of classical d-OPS are given in 
terms of d+ 2 term recurrence relation as well as in terms of functional and differential 
equations. An additional tool in determining the integral representation for such class is 
presented. 

Furthermore, we study the inner products which generate the sequence of orthogonal 
polynomials in the sense of Sobolev. 

Using the theory of Riordan group and d— orthogonal polynomials, we shall show that 
sequences of d— orthogonal polynomials, can be also generated by Riordan group. We 


interpret some families of d— orthogonal polynomials in the frame of Riordan group. 


keyword: 
Orthogonal polynomials , d—Orthogonal polynomials , classical polynomials , Riordan 
Matrix 
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في هده الرسالة قمنا بدراسة كثيرات الحدود د-المتعامدة بمفهوم سوبولاف. وبعد ذكر بعض خصائص 
تعامد سوبولاف e‏ قمنا بإنشاء جداءات داخلية التي تولد متتاليات كثيرات الحدود د المتعامدة بمفهوم 
سوبولاف وأظهرنا أن كثيرات الحدود الكلاسيكية بالنسبة للشعاع دو البعد - د - للتطبيقات الخطية انها 
ليست كلاسيكية بمعنى سوبولاف 


ثم درسنا العديد من خصائص كثيرات الحدود الكلاسيكية د-المتعامدة . أولاء نقدم برهان جديدا على 
معادلة بيرسون باستخدام د- ضعف التعامد و توصيف جديد لكثيرات الحدود الكلاسيكية د-المتعامدة 

باستعمال د- and‏ التعامد. كما تم التعبير عن خصائص أخرى باستعمال د-شبه التعامد و ددضعف 

.التعامد. وبالإضافة إلى ذلك» نقدم أداة لتمثيل التكامل لهذه الفئة من كثيرات الحدود 


والموضوع الثاني في هذه الأطروحة مخصص لكثيرات الحدود د-المتعامدة ومجموعات ريوردان. 
وقد استطعنا أن نبين أن تسلسل كثيرات الحدود د-المتعامدة يمكن أن يتولد أيضا عن هذا النوع من 
المجموعات. وعلاوة على ذلك» عرضنا بعض كثيرات الحدود د-المتعامدة من معاملات هذه 
المصفوفات. النتيجة الرئيسية التي تم الحصول عليها في هذه الدراسة هي أن معكوس مصفوفة ريوردان 
الأسي هو جدول معاملات عائلة من كثيرات الحدود د- المتعامدةء إذا وفقط إذا كانت مصفوفة ستيلجس 
تساوي مصفوفة الإنتاج 


الكلمات المفتاحية 


كثيرات ca gaat)‏ كثيرات الحدود د-المتعامدة» كثيرات الحدود الكلاسيكية»ء مجموعات ريوردان 


Résumé 


Le thème principal de ce travail porte sur les polynomes d— orthogonaux au sens de 
Sobolev. Après avoir cité quelques propriétés de l’orthogonalité de Sobolev, nous avons 
construit des produits Scalaires qui génèrent des suites d— orthogonaux au sens de Sobolev 
et nous avons montré que les suites 2—classiques par rapport au vecteur 7 de dimension 2— 
des formes linéaires ne sont pas definie une suite 2— classique au sens de la généralisation 
de Sobolev 


Nous avons ensuite étudié de nombreuses caractérisations de la d—SPO classique. 
Tout d’abord, nous donnons une nouvelle démonstration de l’équation de Pearson en util- 
isant l’orthogonalité faible ainsi que sur une nouvelle caractérisation du d—SPO classique 
en termes de la d— quasi-orthogonalité. D’autres caractérisations ont également été ex- 
primées en termes de la d—quasi orthogonalité ainsi que de l’orthogonalité faible. En 
outre, nous présentons un outil pour déterminer la représentation intégrale de telle classe 


de polynomes. 


Le second theme dans cette thèse est consacré au polynomes d—orthogonaux et les 
groupes de Riordan. Nous avons pu montrer que les suites de polynomes d—orthogonaux 
peuvent étre également générées par ce type de groupe. De plus, nous avons exhibé 
quelques familles d—orthogonaux à partir des coefficients de ces matrices. Le résultat 
principal obtenu dans cette étude est que l’inverse d’une matrice de Riordan exponen- 
tielle L est un tableau des coefficients d'une famille de polynómes d—orthogonaux si et 
seulement si la matrice Stieltjes est égale à la matrice Production. 

Mots clé: polynómes orthogonaux , d— Orthogonalité , polynómes classique , Matrices 
de Riordan 
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Introduction 


1.1 Orthogonality 


The notion of orthogonal system of functions appeared through the study of certain 
problems of functional analysis (integral equations, strum-Liouville problem and more 
generally, boundary problems in partial differential equations). Orthogonal polynomi- 
als in general, and classical orthogonal polynomials in specific, have been the subject of 
extensive work. They are related to many problems of applied mathematics, theoretical 
physics, chemistry, theory of approximation, and several other mathematical branches. In 
particular, their applications are widely used in theories such as Padé approximants, con- 
tinuous fractions, spectral study of Schrodinger, discrete operators, polynomial solutions 


of second-order differential equations, and others. 


The first family of orthogonal polynomials known as ” Legendre polynomials” appeared 
with the first work of A.M.Legendre on the planetary movements in 1784. Actually, he 
established some common properties and he derived a second order differential equation 


which has as solutions this families. He also studied its zeros. 


Subsequently, other families were introduced. For instance, Hermite (1864), con- 
structed a new family, called Hermite polynomials, which was used in the interpolation 
theory. 

N. H. Abel and V. L. Lagrange P .L. Chebyshev began with T.J Stieltjes the creation 
and development of general theory of orthogonal polynomials by the use of continuous 
fractions. 

E. N .Laguerre introduced in 1897 a family of polynomials (called Laguerre polynomi- 
als) and showed the link between it and the continuous fractions. 

In the middle of the nineteenth century, Jacobi introduced a new family which general- 


ized the polynomials of Legendre and Chebyshev. These three families (Jacobi, Laguerre 
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and Hermite), are known as the "classical polynomials". 

Other authors such as T. J. Stieltjes, L. M. Humbureger, T. Carleman, F. Hausdorff, 
and M .H. Stone have shown that orthogonal polynomials solve also the problem of mo- 
ments via continuous fractions. 

Bochner in 1929, demonstrated that the sturm -Liouville differential equation, 
defined by 


o(z)y" + B(x)y + y(x)y(x) = Ay(x) (1.1) 


has a polynomial solution if and only if deg a(x) € 2, deg f(x) € 1 and (x) = 0. 
Furthermore, by using a linear change of variable, Bochner has also shown that there 
are only 4 families of classical polynomials orthogonal solutions of (1.1) which are Jacobi, 


Bessel, Laguerre and Hermite. 


The second characterization of classical orthogonal polynomials is the Rodrigues 
formula which gives us the explicit expression of the classical orthogonal polynomials in 
terms of a and Û coefficients in equation (1.1). 

In 1938, Krall demonstrated that the classical polynomials satisfied in fact an even 


order differential equation of the type 


2r 
3 ai(z)y = \y. (1.2) 
i=1 


The results of Bochner and Krall have been generalized by Louriero and Maroni in 
2008 and 2010 by giving the link between the coefficients of (1.1) and (1.2). 


Other criteria for the classification of the classical orthogonal polynomials have been 
recently given by Kil. H. Kwon and L. L. Littlejohn (1996) on a Sobolev class, i.e., a 


classification based on a symmetric bilinear form defined by 


b(p,9) = (pq) + (r,pq), p,qa © 2, 


In 2006, Jamei and Koepf gave a class of a differential equation which generalizes the 
Bochner equation and the Sturm-Liouville problem. More precisely, they considered the 
case when the degree of the coefficients œ and Û are > 2 and also when they are rational 
functions. They showed that the solutions of this class are symmetric and orthogonal and 
they gave explicit expressions of their weight function. 

The notion of orthogonality has undergone a great deal of generalization. It began 
with the notion of (1/p), (p > 1) orthogonality (A.Boukhmis 1988), the d-orthogonality 
(P.Maroni 1989), the biotogonality (Brezinski 1992), and finally the multiple orthogonal- 
ity (Aptekarev, van Assche ...). These generalizations intended to make the phenomena 


modeled by differential equations more comprehensible, to improve the approximation 
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methods, and to increase the order of the differential equations having these polynomials 
as solution. 

The notion of d—orthogonality appeared in the thesis of J.Van Isegham [82]. This 
notion was used in the study of the Padé approximation of simultaneous formal d—series. 
Subsequently, it was observed that d—orthogonal polynomials satisfied an recurrence re- 
lation of order (d+ 1). 

After two years, the first properties of a sequence of d— orthogonal polynomials have 
been given in B4]. The author gave a new algebraic approach showing many interesting 
characterization for the d—orthogonality, based on the orthogonality's vectorial form. In 
addition, he introduced the notion of d—quasi orthogonality in the same article. 

Many authors have attempted to improve this theory. Although the classical char- 
acter has been approached by Douak and Maroni in numerous articles, other researchers 
(Y. Ben Cheikh, N. Ben Romdhane, K. Douak, I. Lamiri, A. Bokhemis, E. Zerouki, A. 
Saib, etc.) studied certain problems which lead to the construction of many d-analogues 


of the classical families and to the discovery of new ones. 


In 2006, A. Boukhemis and E. Zerouki found an interesting result where they deter- 
mined a linear differential equation which admits the 2-orthogonal classical polynomials 


as solution given by 
(3) (2) (1) — 
4ر1‎ )2:( Pays (@) + Rna(®) Pays )2:( + Rao(t)Pais(r) + Rna(®)Pata(x) = 0, 
with R,;(z) for à = 1,...,4 are polynomials of degree less than or equal to 7 and 
Ry a(x) = Fra(x)S3(x), 


where S3(x) is a polynomial which does not depend on n and whose degree is less than 3. 
The classical character, i.e. the sequence of orthogonal polynomials in which the 
derived sequence is orthogonal (OPS). This is known as the property of Hahn. 


1.2 Riordan array 


The concept of Riordan arrays has been introduced in 1991 by Shapiro et al with 
the aim of generalizing the concept of Renewal Array defined by Rogers [67], and have 
pointed out its connection to the Umbral calculus [68]. Their basic idea was to define a 
class of infinite lower triangular arrays with properties analogous to those of the Pascal’s 
triangle whose elements are the binomial coefficients. 

The Riordan groups are particularly important in studying combinatorial identities and 
sums. This concept has been investigated by Sprugnoli [75] who pointed out the relevance 


of these matrices from the theoretical and practical point of view. In addition, his work 
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verified that many combinatorial sums can be solved by transforming the generating 
functions. So, we can see that the Riordan array concept is particularly important. For 
example, identities involving Stirling numbers can not be treated by methods related 
to hypergeometric functions. In [76], Sprugnoli paid attention to the identities of Abel 
and Gould. Further applications in reciprocal functions are discussed in [83], on 
subgroups of the Riordan group in Peart and Woan [66], on some characterizations of 
Riordan matrices in Rogers [67] and Merlini et al [60], and in many interesting related 
results in Cheon et al. ZI], He et al [38], Nkwanta [42], Shapiro [74], Barry [7] and so 
forth. 

Successively, some other aspects of the theory have been studied [60]. Further gener- 
alizations of these last characterizations are given in [39] [36]. 

The Sheffer sequence is a very general concept that includes many polynomial sequences 
as special cases. The concept of Sheffer Riordan group of Sheffer-type polynomials is 
defined by He et al. [38]. They prove the isomorphism between the Sheffer group and 
the Riordan group and present an equivalence between the Riordan array pair and the 
generalized Stirling number pair. 

The lower triangular matrices and matrix factorization problems have catalyzed many 
investigations in recent years. Peart and Woodson [64] did some researches on this prob- 
lem. They showed that some classical Riordan arrays have triple factorizations of the 
form L = PCF, where P, C, F are also Riordan arrays. We notice also that the inverse 
of the Stieltjes matrix represents a coefficient of the sequence polynomials [4] [S1]. 
In [79], the authors present new factorizations, and provide extensive examples of families 


of classical polynomials. 


The main results of this thesis 


This thesis contains six chapters. 
Chapter 3 


In Chapter 3, we are interested in the study of orthogonal polynomials of Sobolev. To 
this aim, we first investigate properties of symmetric bilinear forms and in particular the 


one which is defined by 


$(p,q) = (o,p q) + (r.p d), p,q € P. 


We give necessary and sufficient conditions for orthogonal polynomials with respect to 
the symmetric bilinear form © given above, which satisfy the Bochner’s second - order 
differential equation. A classification of the orthogonal polynomial sequences of Sobolev 


solutions of the Bochner differential equation will be investigated also. 
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Finally, we study the regularity of symmetric d—dimensional vector forms © = 
(ġo, Qı, Ris وم‎ defined by 


br(p, q) := (aô(x — c), p d) + (Tr, P q) 


More precisely, we show in the case when, a — 0 and 7 is 2-classical functional 
moment, the sequence of classical 2-orthogonal polynomials satisfies a non-homogeneous 


third-differential equation of Boukhemis type. 
Chapter 4. 


Many of the properties of the classical orthogonal polynomials have nice extensions 
in this d-orthogonality setting: there will be a higher order linear recurrence relation 
[16] B4], there are nice functional equation (Pearson equation) [29] BO], linear combination 
(finite-type relation) between the sequence and its derivatives [58] and differential equation 
[19]. In addition, if the same ideas are to be found, each aspect of the theory is more 
complicated and the links between the different aspects are not so clear. 

This chapter deals with algebraic aspects of classical d-orthogonal polynomials. In 
the first section, we recall some definitions and characterizations of the d—orthogonality 
which we need in the sequel. In addition, it is well known that there is exactly 2° sets 
of d—symmetric d—orthogonal polynomials. We discuss here again this result from the 
combinatorics point of view. The main results are listed in the third section. Many 
characterizations of the classical d-OPS are presented. First, we give a new proof of the 
Pearson equation in the algebraic aspects with the aid of the weak orthogonality as well 
as on a new characterization of the classical d-OPS in terms of d—quasi-orthogonality. 
Further characterizations also expressed in terms of d-quasi-orthogonality as well as on 
the weak orthogonality. Further, we present a set of tools help to determinate the integral 


representation for a such class of polynomials. 
Chapter 5 


In this chapter, we shall show that the sequence of d-orthogonal polynomials, i.e., a 
sequence of polynomial satisfying a linear recurrence relation of d + 2- terms, can be 
also generated by a Riordan array. General references for orthogonal and d—orthogonal 
polynomials can be found in [24 54]. Links between Riordan arrays and orthogonal poly- 
nomials have been explored in [3] Af. Note that there are some examples in the literature 
of Riordan arrays that generate some families of polynomials sequence satisfying a recur- 
rence relation in four terms and more as well as a recurrence with variable coefficients 
[?, [6| B3]. In this chapter, two examples of d-orthogonal polynomials are presented, for 
which an ordinary Riordan array exists. The first one is the d-Chebyshev polynomials 
of second kind. The second one presents in fact the co-recursive d-Chebyshev of second 


kind. The next sections deal with Stieltjes matrix. We generalize the result of Peart and 
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Barry where we show that the corresponding Stieltjes matrix characterizes the sequence 
of d-orthogonal polynomials and we obtain the explicit form of that Riordan array. For 
the exponential Riordan array, we show again that Stieltjes matrix characterizes the se- 
quence of d-orthogonal polynomials. Some examples are discussed. Finally, we improve 


the Sheffer Riordan array by giving an example. 


1.2. Riordan array 


Preliminaries and notations 


The present chapter reminds some concepts about orthogonal polynomials of a real vari- 


able and their fundamental properties. 


2.1 Orthogonality in general 
Let P be the linear space of complex polynomials of one variable and P’ its topological 
dual space. We denote by (u, P) the action of u and PE P. 


2.1.1 Free and monic polynomials 


Let {Pa}, ~o be a sequence of polynomials with coefficients in C. 


n20 


Definition 2.1 The sequence {P,},,9 is said free if and only if deg(P,) = n, Vn > 0 


Definition 2.2 We say that a free sequence {P,},,., is monic if each polynomial P, is 


written in the form : 
n—1 
P,(z) = a" + 3 dua”, for n > 0 
k=0 


Proposition 2.1 If {P.},>o is monic, then there exist a unique sequence {B,} and a 
unique array Xnv, 0 X v € n, such that 

F(z) = 1, Pix) = v — Bo 

P,+2(2) = (x ME 8) Pa+1 m Pt oco y) = 0 
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Definition 2.3 Let {P,},:, be a sequence of free polynomials, the function G(x; t) which 
can be developed into a power series of 1 is called the generating function of the sequence 


{Pr} nso if it is put in the form : 


Gat) = S Paat where c, عو‎ 0 ,Vn > 0. 


n20 


2.1.2 Regular and definite positive forms 


Definition 2.4 We call the dual sequence of the free sequence {Pp} the sequence of 


n20 


linear forms (u,) defined by 
un( Pm) = (un, Pn) = noni n,m > 0, 


where (.,.) is the duality hook. The linear form wo is called the canonical form of the 


sequence {Pa}, sg. 
Proposition 2.2 i) [Z/|For each monic sequence, the dual sequence exists and is unique. 


ii) In this case, we have 


Pn — Unt Pn (2)), n > 0, 


(2.1) 
Xnv = (Uu, Papi (x)) ; 0 € v € n. 


Definition 2.5 The form £ is called regular (or admissible, or quasi-definite) if there 


exists a free sequence {P,},:, such that 


£(P Pn) = 0, 8 Em, 


(2.2) 
LP) عو‎ 0, n>0. 


Such a sequence is called orthogonal (or regularly orthogonal) relative to £. Such a 


sequence is free and unique up to a multiplicative constant. 


Lemma 2.1 For any £ € P’ and any integer p > 1, the following statements are equiv- 


alent [54] 55: 
I) (£, Pp-1) # 0, (és Pn) = 0, n > p, 
ID 20, EC, 0<v <p-1, رد‎ £ 0 such that 


2-1 
£=) has 
Lemma 2.2 The sequence {Pa}n>o is orthogonal with respect to £. Necessarily, £ = 


A£o, A #0. In this case, we have [54 


= (A PA) P uo, n> 0. 


when £ is regular, let A be a polynomial such that A£ = 0. Then A = 0. 
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Remark 2.1 For any linear form £ and for any polynomial 7, we can define the two 
forms D£ = £' and r£ by 


(Lf) = (Lf), ef) =(£,0f), TEP. 
and thus, the dual sequence {ŭn} >o of {Qn}, is given by 
à, = - (n + 1( ريبس‎ n2 0, (2.3) 
where Qn (x) =(n+1) Pla (2), n > 0. 
Definition 2.6 The sequence {£,},,-, defined by 
£n := L(x"), n >0 


is called a sequence of moments of the form £. We say that the form p is real if the 


moments £n, n € N, are real. 


Theorem 2.1 49] The form £ is regular if and only if the Hankel determinants 
defined by 
An = det(£;,;)7; 9, 7 > 0 


are all non null. 


Proposition 2.3 The following conditions are all equivalent 
1) {Pa}nso ds orthogonal polynomials sequence relative to p. 
2) 1 

0 if degr «n, ع"‎ <1 


Z 0 if degrî =n‏ = ((ند)رظ(ته)ع) ده 


3) 
Zum PE) = AnOnm, An 4 0,0 € men 


Definition 2.7 A form £ is said positive definite if (£,7(x)) > 0 for any polynomial 
m = 0 such that a(x) > 0 for all real x. 


Proposition 2.4 If a form £ is positive definite, then it is regular and all mo- 
ments of £ are real. 


Theorem 2.2 [49] A form £ is positive definite if and only if it is real and the 
determinants of Hankel A, are strictly positive for all n > 0. 
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2.1.3 Fundamental recurrence. 


Theorem 2.3 [49/Let £ be an admissible linear form on P and {Pn},>o the orthogo- 
nal monic sequence relative to £. Then this sequence satisfies a linear recurrence of order 


two, i.e. there exist two series of numbers {An}n>0 and {Bn}n50 such that 


(2.4) 


P(x) = 1, Pi(z) = x — Bo 
Past) = (£ — bn) Pea = Ona a) , n>0 


with the regularity condition a, # 0 for all n > 0. Moreover, if we set P,(x) = x" + 


np de 
b," +... for all n > 1, we have, Bn = 6,43 — bn+1 and an = where A, =1 


(by convention). 


Theorem 2.4 Let {Pr}nso be a monic sequence of polynomials satisfying the recur- 
rence [2.4] where an #0 for any n > 0.Then there exists a single linear form £ such that 
£(Po) = ao and £(Pm(£)Pa(£)) = 0 jor ang n عو‎ m. 


Corollary 2.1 Let {Patn>o be a monic sequence of polynomials satisfying the recurrence 
Then the form u is positive definite if and only if 5, € R anda, > 0 for any n > 0. 


Definition 2.8 Let {P,}, 
rence We call associated polynomial sequence of {Pa}n>0 the sequence denoted by 
UP, c) aso, for any c € R, defined by : 


be a monic sequence of polynomials satisfying the recur- 


Po(z,c) = 1, Pis; 6) =g — Be 


Pa+2(2,€) = (£ — Bateti)Pnti(Z,c) = Onte+1Pn(Z,e) , n >0 


Theorem 2.5 Let {P,},59 be a sequence of orthogonal polynomials with respect to 
the positive definite form u, then each polynomial P, admits n real and simple roots for 


alln > 1. Moreover, the zeros of P, and P,:1 are alternate, i.e. 
Ln+1,v < nw < Ln+1,0+1 < YXny41; N > 1,v > 1, 
+1 ; 
where (14,,),. ., and (x41, v); 4 are respectively the zeros of P, and Pa41. 


2.1.4 Classical sequences 


Definition 2.9 Let {Patn>o be a sequence of orthogonal polynomials. This sequence is 


called classical if these sequences of successive derivatives are also orthogonal (Hahn's 


property). 
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Theorem 2.6 The classical sequences are the only real polynomial sequences 
solution of the second-order differential equation of the hypergeometric type of Sturm - 
Liouville 

a (2) y" + B (z) y + ny = 0, (2.5) 
where a and B are polynomials of degree two and one respectively and where A, is the 


following constant: 
An = o" (zx)n(n — 1) + nf'(x) 


Theorem 2.7 The classical orthogonal polynomials are given by the formula of Ro- 
driguez 
کے = ميم‎ ۳ (a" (2)u(r)) (2.6) 
nib] ss e , . 
knw(x) dz” 


where kn is a non-zero constant and w is a solution of the differential equation (aw)! = pw. 


Corollary 2.2 If £ admits the following integral representation 


b 
£(P)= | p(z)w(x)dx , for any polynomial P, 


then the sequences of polynomials given by [2.0 are orthogonal to the weight function w. 


2.1.5 Integrated Functional Representation 


We know, from Duran’s classical theorem concerning the problem of moments of Stieltjes, 
that every moment function can be represented by a measure in the form of an integral 
of Rieman-Stieltjes. Now it is interesting to establish the integral representations of the 


moment functionnelle. 


Definition 2.10 Let @ be a function defined on J C R. 6 is said to be a rapidly 


decreasing function on / if, for every pair of integers (n,m), we have 


m 


iir (zr) > مهل‎ for any vel. (2.7) 
qm 


Definition 2.11 A bounded non -decreasing © whose moments 


wi Û rao) 


oo 


are all finite is called a distribution function, 


Theorem 2.8 [34/ For any sequence {Un},>o of complex numbers there exists a fast decay 


function on Rt (respectively on R) such that 


+00 us 
Un E x"ó(x)dx respectively ju, HI x'p(x)dz. (2.8) 
0 3 


OO 
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Theorem 2.9 [/6) Let {Patn>o be a regular orthogonal polynomial sequence satisfying 
the differential equation [2.5] If w(x) is an orthogonalization weight function of {Pa}n>o 


then w(x) satisfies the following differential equation 
(a(x)w(z)) — B(x)w(x) = g(x) (2.9) 
where g(x) is a weight distribution with all its zero moments, i.e. 
(glx), z”) =0,n > 0. (2.10) 
Conversely, if w(x) is a distribution satisfying the following conditions 


I) w(x) decreases rapidly when |x| tends to oo; then (w(x), x") exists and is finite for all 
n 20; 


II) w(x) defines a non trivial moment function, i.e., its moments are not all null ; 


III) w(x) is a solution of equation in the distribution space in R ; 
then w(x) is said orthogonalization distribution of {Pn}, >o: 


Definition 2.12 Relation is called the non-homogeneous weight equation for the 
differential equation 


2.2 The d—Orthogonality 


Let us consider linear forms To, ..., la 1 (d > 1). 


Definition 2.13 The sequence {Pabnso is called d-orthogonal polynomials se- 
quence, in short à d-OPS, with respect to the d—dimensional vector of linear forms 
T = (Ts, ..., Da i)?, if it fulfills 

Tas” Pa (xz) =0, n> md+a+1, m >0, 


(2.11) 
)1' ”تتبن‎ Pee #0, m >0, 


for each 0 € a € d — 1. 


In this case, the d-dimensional functional — is called regular. 
Remark 2.2 [18]1 If (P), is a d—orthogonal polynomial sequence, then its polynomials 


are exactly of degree n and can hence be normalized , thus it follows the uniqueness of 
the sequence. 

2) We say that the functional T is regular of dimension d 

3) When d — 1, we find the usual notion of regular orthogonality 
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Theorem 2.10 [54/It is necessary and sufficient, for the functional T of dimension d to 
be regular, that the following determinants be non-zero : 


lo ri ses Péri 
I I> as Pindto 
Hma+o = [li I diu née , m20, v=), 0-1 (2.12) 
Dh. ls ds vane ee 
Vin Tm ee haat 


We have posed Fp = T(z”) = (T(a"), T(z"), ..., P*(z7))7, n > 0. Represents actually 
d line. 


Therefore, each determinant Hma+v is of order md + v if md+v > 1. 
By convention Ho = 1. Each of the first m lines of the determinant [2.12] represents 
in fact d line. 
we have (25) : 


m-—1 d vU 
Hmaty = I] I] D Phare) I] Ue Pmara-1(7)) 
u=0 a=1 a=] 


for any m > 1 andl<v<d-1 


m—1 d v 


Hma = [[ [[U*(" Ree) ||] 21 


and 
aln Hasta 
PP nat للد‎ n> 0,1 > به‎ > 0 
Let us now recall some characterizations we need in the sequel. 


Theorem 2.11 Let {Pr}nso be a monic sequence of polynomials, then the following 
statements are equivalent. 


(a) The sequence {P,},,50 is d-OPS with respect to T = (To, ..., Dos. 
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(b) The sequence {P,},,9 satisfies a recurrence relation of order d +1 (d > 1): 
Pm+a+1 (x) = (x — Ûm+a) P, He Se yl m+d—1—v (x), m = 0, (2.13) 
with the initial data 
ر(ته) مر‎ =1, P(x) =z — bo, (2.14) 


Pa) = (a= Baa) Poa) د‎ So Ps un), 2<m<d 
and the regularity conditions رج‎ #0, m > 0. 


(c) For each (n,v),n > 0,0 <v € d — 1, there exist d polynomials Qf ,, 0 € u <d—1 


such that 
0-1 
gS bh ty, n>0, USE )2.15( 
p=0 
and verifying 
deg D, = n, 0<v<d-l, and if d 2, 
deg dit, <n, ü stus pe if 1<v<d-1, (2.16) 


-0 > بع 0 if‏ ,1 -0 > بر عه 1 د degon,Sn-l,‏ 


Notice that the coefficients {8,,}and {Ym4i}are given by [54] 


Da = (sP), 0 > نه‎ > 0-1 (2.17) 
Me = (Untu, Para), OSU <S 0 - س1‎ < 0 


Proposition 2.5 For any d— OPS satisfying the recurrence relations [2.13] and [2. 14] 
the first moments 

(uo)i — fo 

up = 5ق‎ + y 
2 = Po Pi 
3 = Bot 61 + Ba 


= 


1 


e 


) 
) 
2)3 = 
u3)4 = E + 61 + b2 + Bs 
1)a = 67 + B2 + B + bobi + +f + 
و(‎ = 8 + f (260 + 81) + YE? 
) 
a= 
a= 


LOS NNN SSS 


uo 

U2)4 = AH ARE Ra ce LH + کہ‎ 

uila = e H + qf Bo + 281 + £2) +70? + 
ua = 84 + [265 + (Bo + £1)? + yet + Y "| + 1 (280 + Bı Bo) + 407? 
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Proposition 2.6 [56] The statements in theorem are equivalent to the following 
(d) Xn+d—1,v = 0, TL > U + 1 and Xv+d—1,v A 0, U > 0, 


(e) 


d—1 
LUn = Uni + nun + b» Xn+unUn+i+v, n 20 (ui = 0) 

v=0 (2.18) 
Xn+d-1n Æ 0, n = 0. 


Corollary 2.3 [72/The multiplicity of zeros of any d — OPS is at most d. Moreover, 
any d + 1 consecutive polynomials as well as an d + 1 consecutive polynomials of 
r—associated sequence {Pn}, 59 have no common zero. And any r > 0 the polynomi- 
als PP, PP, p(r*à) have no common zero 


2.3 Quasi - monomiality 


We give some results of quasi - monomiality in order to get the dual sequence of a given 
polynomial set by using the lowering operator associated with the involved polynomials. 
The results obtained will be applied to Boas -Buck polynomial sets. For an integer j, we 
denote by AU? the space of operator acting on analytic function that augment (resp-reduce 
) the degree of every polynomial by exactly j if 7 € N (resp j € Z`), Z being the set of 
non positive integers. This tacitly includes the fact that, if o € AW), then o(1) = 0. 


Definition 2.14 Let c € ACD. A polynomial set {Pn},>o is called a sequence of 
basic polynomials for ø if : 


(i) Pela) =1, 
(ii) P,(0) = 0 whenever n > 0, 
(iil) oF, (2) = nP a). 

As a consequence of this definition, we mention the orthogonality relation 


o” Pal0) = n!őnm, “HS ou. (2.19) 


Starting from a polynomial set, We give a method of constructing its dual sequence in 
terms of its lowering operator 


Lemma 2.3 [8| Every o € AW has a sequence of basic polynomials 
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Lemma 2.4 Let {Qn}nso be a polynomial set and let © be its lowering operator. Let 


{Pr} nso be a sequence of basic polynomials for o. Then there exists a power series 


p(t) = oa Qo + 0, 
k=0 
such that 


gie = Pn, de 


Let us combine the two preceding lemmas, gives us the following theorem 


Theorem 2.12 [8/Let {Qr}nso be a polynomial set and let © lowering operator , and let 


{£n}n>o be its dual sequence. Then there exists a power series 


e(t) = 9 ant”, os z 0, 
k=0 


such that 


(Ln f) = Slo" Neo = TESI n=O feP  — (Qm) 


n! 
If we apply Theorem we obtain a further method to construct the sequence 
{£n وجو[‎ for general case from which we deduce the following expansion theorem 


Theorem 2.13 Let {Qn},>g be a polynomial set and let o be its lowering operator , 
and let {£n} >o be its dual sequence. Then there exists a power series 


y(t) = 3 ant”, 20 > 0 
k=0 


such that every analytic function f has the expansion 


oo 


f) = 3 C90), (y (2.21) 


If, morever, the translation operator Ta commutates with o, then 


Jeta = Y 72000, (2 


Corollary 2.4 Let {Qn}nso be a Boas -Buck polynomial set generated by the formal 


relation 


deg- dou » Onl) in (2.22) 
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where 
A(t) = Sat", B= bit* and C(t) = V o? (2.23) 
k=0 k=0 k=0 


are three formal power series with the condition aocobx # 0, for all k. Let v := v, € 
AC) such that 


vB(xt) = tB(xt) 


Put o = C*(v) where C*is the inverse of C i.e. 


C* (CŒ) eO) = t, with C* (t) = Y arm, c5 #0, 
k=0 


then every analytic function has the expansion 


n! 


= 5 PLUS 


If, moreover, the translation operator T, commutates with o, then 


n! 


fG@+a)= 5 E 


Definition 2.15 [9] A polynomial set {P,},,:, is called quasi -monomial if and only if it 


possible to define two operators c and 7, independent of n, such 


cP (fo) = Pa) (2.24) 
TP) = P(x) 


Theorem 2.14 [Y Every polynomial set is quasi-monomial 


Theorem 2.15 [9| Let (o,r) € AC? x AM and let P :— {P,},>o be a polynomial set 
generated by 


G(z,t) = 3 Pala) yn (2.25) 


Then we have the equivalences: 


oG(z,t) = tG(x, t) جح‎ oP, = nPa1 


and 
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باط حم نر حم TG(z,t) = ar)‏ 


where o :— o, and T := Tz 
Definition 2.16 |?] the operator ø is called lowering operator of {Pn(Z)},>o 


Remark 2.3 According to Theorems and , we deduce that if a function of 
two variables G(x,t) has an expansion of type when deg P, = n, there exist two 
operators (o,r) € AT) x AU) such that G{(x,t) is a eigenfunction of o :— o,,associated 
with the eigenvalue t and is also a solution of that equation associated with the operator 


T. 


2.4 Riordan arrays and orthogonal polynomials 


Riordan arrays give us an intuitive method of solving combinatorial problems, and help 
to build an understanding of many number patterns. They provide an effective method 
of proving combinatorial identities and solving numerical puzzles as in [61] rather than 


using computer based approaches rq 
In this section we review mainly known results related to integer sequences and Riordan 


arrays that will be referred to in the rest of the work. 


2.4.1 Integer sequences and generating functions 


Definition 2.17 The ordinary generating function (0.g.f.) of a sequence a, is the formal 


power series 
g(x) = + anz” (2.26) 
n=0 
Definition 2.18 The exponential generating function (e.g.f.) of a sequence a, is the 
formal power series 


h(x) = 3 an (2.27) 


Definition 2.19 The bivariate generating functions (b.g.f./s), either ordinary or expo- 


nential of an array (a@,,,) are the formal power series in two variables defined by 
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alz, y) = 25 An x" y" (o,g,f) and (2.28) 


n,k 
x” 
b(x, y) = p» ane TU (e,g,f) 
n,k 1 


Notation 2.1 The coefficient of x" is denoted by |x"]g(x) and from the definition of the 
e.g.f., we have n!|x"]h(x) = (Rx) We refer to the inverse of f as the series reversion 
: n! 


f: 


Lagrange inversion [85] provides a simple method to calculate the coefficients of the 


series reversion. 


Theorem 2.16 Let olu) = oku" be a power series of C[[u]] with do # 0. 
Then, the equation y = zó(y) admits a unique solution in C[[u]] whose coefficients are 
given by 


y(2) = قم‎ Un = 
n=0 


The Lagrange Inversion Theorem may be written as 


"|G (Fe) == [7] Ga) £2] 


The simplest case is that of G(x) = x, in which we get 


pee L fa] ey 


If the product of two power series f and g is 1 then f and g are termed reciprocal 
sequences and satisfy the following. For o.g.f./s we have [161] 


Definition 2.20 A reciprocal series g(x) = X ga,z" with مه‎ = 1, of a series f(x) = 
ورور ا‎ bng” with bo = 1, is a power series where g(x) f(x) = 1, which can be calculated as 
follows 


3 anz” = — 3 »» و0 ا‎ ag = 1 (2.29) 
n=0 


n=0 i=1 
and for e.g.f.’s we have 


Definition 2.21 A reciprocal series g(x) = So di with ag = 1, ofa series f(x) = 
n! 
D ba with bọ = 1, is a power series where g(x) f(x) = 1, which can be calculated 5 
n! 


follows 
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oo gh co n x” 
MEET = 2 D bey; ag = 1 (2.30) 
n=0 n=0 i= 


2.4.2 The Riordan group 


Definition 2.22 A Riordan array is a pair (g(x), f(x)) in which g(x) and f(x) are formal 
power series such that g(0) # 0 and f(0) =0. 
The pair defines an infinite, lower triangular array L = (l;, دعي مزلم‎ where 


Ine = fr] g(x) f£ (2))* 
If f'(0) Z 0 the Riordan array is called proper. 


Remark 2.4 From this definition, it easily follows that g(a) f(x)" is the generating func- 


tion of column k in array 


Remark 2.5 A lower triangular matrix L is a Riordan Array if the generating function 
kt” column is 


a Fi Tor os 01:93. 
with 


where the generating function of the two column vectors are A(z) = J} anz” and B(z) = 
n20 

X bz" respectively. The identity is true if and only if the following equation holds : 

n20 


g(z)A(f(z) = B(z) 
Proof. We look at the Riordan Array (g(z), f(z)) column by column, and multiply it by 
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the column vector on the left-hand side 


ao 
ai 
9 gf gf . .| | a 
This yields 
aog + agf + aag f^ +... = glao + af + af? +...) = g(z).A(f(2)) = B(z) 


and we have our result m 


Definition 2.23 The Riordan group 


pa) OC) fl) / GG), f(2)) is a Riordan array and f(x) = fo + fiz + faz? +, 
7 where fo = 0, fı = 1 


l.e., 
each member of R is a lower triangular matrix with 1’s on the main diagonal.The 
multiplication in R is 


The identity is I = (1,2). 
The inverse of (g(z), f(z)) is ( 
of f(z), i.e. 


a 
A 
up 
z 
a 
D 
ب‎ 
6 
E 
G 
= 
ct 
E 
© 
Q 
© 
z 
md 
© 
D 
et. 
© 
=) 
ES 
3 
< 
© 
ب‎ 
gi 
6 


To check the inverse property, we compute 


i. = 1 0 F 
(zx Fle) o rn = (O) = (1,2) 


2.4.3 The A—sequence 


The A—sequence introduced by Rogers (1978) characterizes the column elements of after 


the first column 
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Theorem 2.18 An infinite lower array D = (dn,k)n,keN, is a riordan array if and only if 


a sequence À = {an} existe such that for every n, k € No it is true 


n€No 


n—k 
dn+1,k+1 = 1 aidn kti 
i=0 


even more, if D = (d(x), h(x)), then the generating function of the A— sequence is such 
that A(x) = D, ax satisfies the equation 


h(x) = zA(h(x)) = A(x) = مك‎ 


2.4.4 The Z—sequence 


The Z—sequence introduced by Merlini et al characterizes the elements of the first column 


of a proper Riordan as follows 


Theorem 2.19 Let D = (d(x), h(x)) = (dyx)nxen. Then a unique sequence Z = 
((Zo, 21, 22, ...( can be determined such that every element in column 0 excluding the ele- 
ment in the first row can be expressed as a linear combination of all the elementin the 
preceeding row with the coefficients identified as the elements of the sequence Z satisfying 


the relation 


n 
dn+1,0 = b3 Zidni n € No 
i=0 
the generating function of the Z—sequence satisfies the equation 


dy NEMORE 
سيد‎ SUA) P Ros ( xi) 


2.4.5 Stieltjes matrix 


In the context of Riordan arrays, We see the Stieltjes expansion Theorem in [65], defined 
as follows. 


Definition 2.24 Let L=(l,x)1# be a lower triangular matrix with l;; = 1 for all i > 0. 
The Stieltjes matrix Sy associated with L is given by Sy, = L-!L where L is obtained 
from L by deleting the first row of L, that is, the element in the n” row and k** column 


of L is given by Ing = ln+1,k 


Using the definition of the Stieltjes matrix above [65] leads to the following theorem 


relating the Riordan matrix to a Hankel matrix with a particular decomposition 
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Theorem 2.20 [65] Let H = {hak = an+k} be the Hankel matrix generated by the 


n,k>0 
sequence 1, a1, ,و0 ,و0‎ ... . Assume that H = LDU where 
1 0 0 0 0 
ho 1 0 0 0 
Lo loi 1 0 0 
L= l mE 2 : 
Un) o l30 131 l2 1 0 
lag lai la2 l3 1 
do 0 0 0 
0 d 0 0 
D=|0 0 d 0 -:-| , ,ل‎ #0, U= LT 


Then the Stieltjes matrix Sr is tridiagonal, with the form 


Bo 1 0 0 
a, & 1 0 
0 Q2 Bo 1 
0 0 ag fs 
0 0 0 «a 
where 
_ dky 
fo = a1, 01 = di, Bk = lxyir — lk,k+1; Oki = J’ k>0 
k 


Now, we state other relevant results from this section, relating to generating functions 


which satisfy particular Stieltjes matrices. The first result relates to o.g.f.s 


Theorem 2.21 [63] Let H be the Hankel matrix generating by the sequence 1, a1, a2, aa, … 
and let H = LDL" .Then Sr has form 


a 1 0 0 
o1 DB 1 0 
0 a 6 1 
0.00 6 
0 0 0a 
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if and only if the o.g.f. g(x) of the sequence 1, ai, a2, a3, .. is given by 


(x) : 
a = 
9 1 — az — œz f 


where 


f(x) =a(1+ Bf + of?) 
Solving both equations above give us the required result. Similarly for e.g.f.s we have 
the following result. 


Theorem 2.22 [65/Let H be the Hankel matrix generating by the sequence 1, a4, az, aa, ... 
and let H = LDL" .Then Sr, has form 


Bo 1 0 0 
a, B1 1 0 
0 ag B9 1 
0 0 os fs 


0 0 0 ag 


if and only if the o.g.f. g(x) of the sequence 1, a1, a2, a3, ..is given by 


qx) = f(a - يفريه‎ with g(0) =1 
where f is the solution of 


% - 1+ 6f tof", f(0) =0 


The proof again in involves looking at the form of the n” column of the Riordan 
array. However, intuitively this result can be from looking at the form of the matrix 


equation L = LS. In the case that L = [g(x), f (z)] is an exponential Riordan array, we 


have the following 
d 


Proposition 2.7 /?] L = ul 
x 


Proof. 


d eo x” 0e gr} 
£x») = 2 973 


n=0 n=1 


OO gr 
= , nci 
n! 

n=0 
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Equating the first column of matrix Land LS we have 


(gr) = Bog(x) + o1g(x) f(x) 


and second column equate to 


L (Fle) = AG) + agla) f Gn 


which gives the required result. m 

The Stieltjes matrix as we have seen above is a tridiagonal infinite matrix which is 
associated with orthogonal polynomials, and therefore have a generalization of the Stieltjes 
matrix to the Riordan group. Referred to as a production matrix jt is defined in the 
following terms 

Let P be an infinite matrix (most often it will have integer entries). Letting ro be 
through throw vector rg = (1,0,0,0,....), we define r; = r; 4P, where i > 1. 

Stacking these leads to another infinite matrix which we denote by Ap .Then P is said 
production matrix for Ap . If we let UT = (1,0,0,0, ...) then we have 


and DAp = A,P where D = (6%;,;41)ij>0 in , P is called the Stieltjes matrix 


associated to Ap 


2.4.6 the production matrix 


In the context of ordinary Riordan arrays, the production matrix associated to a proper 


Riordan array takes on special form: 


Proposition 2.8 [26] Let P be an infinite production matrix and let A, be the matrix 
induced by P. Then Ap is an (ordinary ) Riordan matriz if and only if P is of the form 
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39 0 0 0 مه مي 
a, a 0 0 0‏ & 
a, ao 0 0‏ وه £0 
A 0‏ 04 وبه P= |& o3‏ 
O3 Q2 a, Q‏ يه f4‏ 


Moreover, columns 0 and 1 of the matrix P are the Z— and A— sequences, respectively 
of the Riordan array Ap 
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Sequences of the orthogonal polynomials of 
Sobolev 


3.1 Definitions and properties 


We are given 2 functional moments c and 7 and the bilinear form @ defined by 


o(p, 4) = (o,p رو‎ + (n.p d), p.a EP. (3.1) 
Here, we prove that, when Q is regular, we can classify the polynomials solution of the 
Bochner’s differential equation. We point out that ¢ is symmetric ((p,q) = é(q,p)). 
Definition 3.1 [MI] The numbers 
Onm = Q(x”, 1"), m,n > 0. 
are called the moments of ف‎ (.,.) , for any symmetric bilinear form ¢(.,.). 
Definition 3.2 [ZI] A symmetric bilinear form $ is said to be quasi-definite (resp 
defined positive) if it satisfies the Hamburger condition 
An(d) =: det [dij]; -0 Z 0 (resp. An(%) > 0, n > 0). (3.2) 


Lemma 3.1 Let @ be a symmetric bilinear form. For any normalized sequence of 


polynomials (V, (x)),., we have 


A, (6) = det {6(W;, Y) o ,7 > 0. (3.3) 


ij-0 ? 
Proof. For any n > 0, let A= (84); i-o be the square matrix of order (n + 1) and the 
sequence [9,6] defined by 
j-0 
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Then, we have 


So, 


det[ó(9, &)| _ = (det A)? det [93,9 


i,j=0 ` 


4,7=0 
If we choose a lower triangular matrix such that aj; = 1 and take V (ar) = gi, for 
j = 0,1,2,... n; then 
An(¢) = det lola, e := det [64,5]; io — det [ó(V ;, NE, esa 
Hence, [3.3] holds. m 
Lemma 3.2 A symmetric bilinear form is quasi-defined (resp defined positive) if and 


only if there exist a sequence of polynomials {P,(x)},>-0 and real kn ¥ 0 (resp kn > 0) 
for n > 0 such that 


(Pn, Pn) = knônm n,m > 0. (3.4) 


In this case, an SP {P,(x)},=0 satisfying [3.4] is uniquely (to a non-zero multiplicative 
factor) determined. 


Corollary 3.1 The following statements are equivalent. 
a){Pn(7)},>o ûs a sequence of orthogonal polynomials with respect to @ 


olz”, Pt) = 0 J0em:csSm-.l, 
x”, Pa(x)) # 0. 


b) For every polynomial n(x) of degree m, 


ln(z), Pa(x)) = 0 if 0 > 1-س > ور‎ 
ln(x), Pal) £0 if m=n. 


We now investigate the regularity of the form @ defined by 3.1] 


Remark 3.1 The integral representation of a quasi-definite moment function is always 
possible according to Duran’s classical theorem. It is convenient to define the orthogo- 


nality following in the sense of Sobolev. 
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Definition 3.3 We say that (P,(r)],.9 is a sequence of orthogonal polynomials of 
Sobolev (SOPS) if it is orthogonal to Q defined by 8.1] i.e., 


(Pay Pm) = (o, Pa Pn) + (7, PLPL Y = knônm, with kn #0, n,m > 0. (3.5) 


Lemma 3.3 Let ¢ be a quasi-definite bilinear form defined by [3.1 and let 
{P,(x)} o be an SPOS with respect to ¢. Then, o is the canonical functional of 
(Pa) oies 

(0,1) عو‎ 0.010 (o, Pn(z)) =0, n>1. 


We give necessary and sufficient conditions for the form © to be regular. To this end, 
we investigate three cases. 

Case 1: o is regular. 

Case 2: ع‎ is regular. 

Case 3: o and 7 are not regular. 


Definition 3.4 The sequence of derivatives of the sequence of polynomials {Kn (x, y)}n>0 


is defined by 
d'*$ 


2 _K (ey), n>0. 
pup ee m 

Let us denote by{Qn(x)}, {Ra(x)} 0 and (P,(x)) 9 the corresponding SPON veri- 
fying the orthogonality conditions 


(T, Qu(z)Qu(r)) = kn (T) bmn n,m > 0 
(d, Rm (2) Ry (Z)) =k, (9) Ôm,n n,m > 0 


CLADE = kn (0) m,n n,m 2 0. 


Now, we begin with the first case (c is regular). where ( ô., p ) := p(c), ce R 


Lemma 3.4 [49] The bilinear and symmetric form @ defined by 
(p. 9) := (o, pq) + a(d.,p'q) p,q € P 
is quasi-definite if and only if 


1raKÜP(ec) Z0, n>1. (3.6) 


aD, (z) c 9 (7, €) 
1+ aK (e, c) 


for n>0. (3.7) 
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Moreover, with K_1(x,y) = 0 


or 


R! (x) = 20) n>0 
Fe) Vn > 0 


P(x) = Qo(x) et Pi(x) = Qi(z). 


Now, we consider the second case (c = ad, ). 


Lemma 3.5 [/9| Let ó be the bilinear and symmetric form defined by 
(p,q) = (be, pq) + (7, pq). 


Then, Q is quasi-defined (positive definite) if and only if a # 0 and T is quasi-defined 
(resp. a > 0 and T is positive definite). 

In the case where ġ is quasi-defined, we denote by(Q, (x)) and (Rs (x)) y respectively 
the corresponding SPON verifying the orthogonality 


(T, Qm(t)Qn(x)) = ks (T) Öm,n 


(9, Fo (x) o ()) = kn (0) 9s 


then 


R,(c) =0, n>1 where ô,(x) is a canonical functional of {R (£) 


n=0 
kaa(0) =(n+1)kn(r), n <0 


S = ae 


n +1 n+1 


Ba) = ةل‎ nQn-1(z)dz, n>1, Roz) =1 


Proposition 3.1 Suppose that a > 0 andr is a moment classical functional verifying 
Pearson functional equation 


((x)r) = h(x)r, 
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or 
L(x) = us, i= 1,2. 
1-0 
Then {R,(x)}satisfies the non-homogeneous equation of order two 
In(a)y" + (h(x) — l2(z))y' + (nl3 (z) — Any(@)) = vn; (3.8) 


with 
An = loon(n = 1) -+ lun 


Un —la(c)R,(c) + (Lle) + hle) رط‎ (c) form > 0. 


3 Case c, 7 are not quasi-definite and © quasi-definite. 


Example 3.1 [49] Consider the 2 functional moments c and 7, defined by 


1 1 
T= 0 = gos); 


1 
doo - 3 =] et 0y كت‎ =, >1, 
and oo = 3, oi ا‎ En e 
or (05,9 is the sequence of moments of FM o. 
As 
| 0 
— n= 
n 8’ 
T (z^) = 
0, n>1, 


we deduce that 7 is not quasi-definite. 
In addition, we have A, (c) = 0, therefore o is not quasi-definite. 
Let us now define the bilinear form ¢ by 


1 LENA 
(p,q) := (o, p q) + 3o, P 4), pP, qEP. 


We know that ; 


$(.) = | مهم"‎ + (0 > o 
0 
then, Q is positive definite > 6 is quasi-defined. 


We can construct more such examples via 


Lemma 3.6 [49] For any symmetric bilinear form W, there is a positive-definite moment 


functional o such that the symmetric bilinear form © defined by 


ó(p.q) := V(p', d) + (o, pa) 
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is positive-definite. 
In particular, for any FM T, there exists a positive- defined FM o such that 


o(p, 9) := (o,p q) + (T, pq) 


is positive-definite. 


3.2 Sobolev Orthogonal polynomials and Second or- 
der Differential Equations 


Remark 3.2 Let م‎ be a quasi-definite bilinear form given by the relation where > 


and 7 are two FM. 
Our aim here is to determine necessary conditions for the SOPS with respect to 6, 


satisfy a second order differential equation having the form 
Ly] (x) = la(zx)y" + h(x)y = Any, with L[] the differential operator (3.9) 


where the coefficients ,نا‎ i = 1, 2 and An, n > 0 are as follow 


An = n(n = 1)l22 + nl, n > 0 


l? + l2 Æ 0. 


Theorem 3.1 Consider the symmetric bilinear form © defined by [3.1] The fol- 
lowing statements are equivalent. 


a) The differential operator L |.|, defined by [3.9] satisfies the following relation 
olL [v], a) = o(p, Lidl), p,q € P. (3.10) 
b) The two functional moments o and v satisfy the following functional equations 
(la(z)o)' —l(x)o = 0 (3.11) 


lo(z)r' — h(zx)r = 0 (3.12) 


c) The moments of ¢, o and T are given by 


Omn = Omin + MNTm+n-2, M, N >0 


(3.13) 


Tog — Tas 0, 
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(nlog + بره( نا‎ + (nla + ho)on + nlages 1 = 0, n > 0, (3.14) 
04 = 0 | 
and 
| [(n + 2)23 + hi] ua + [(n + ler + hol T + nlooTn-1 = 0, n > 0 (3.15) 
T—1 = 0 | | 


Moreover, if ó(.,.) is quasi-definite and (P,(z)), y isa SOPS, then the conditions 
a), b) and c) are equivalent to condition d), given below 
d) {Pa(x)} o satisfies the differential equation [3.9 


Definition 3.5 The functional equations and are the Sobolev weight equations 
of the operator L|. with respect to the form ¢(.,.). Equations and [3.15] are the 
moments of Sobolev for L |.] with respect to ¢(.,.). 


Remark 3.3 Theorem B.1]is a generalization of the Krall theorem concerning the differ- 


ential equations of the 24 order of the series of orthogonal polynomials. 


Corollary 3.2 Let $(.,.) be a symmetric bilinear form defined by[3.1] Then there exists a 
sequence of orthogonal polynomials of Sobolev with respect to $(.,.) satisfying a differential 
equation of type [3.9] if and only if both of the following are satisfied 

1) $(.,.) is quasi-definite 

2) The moments {os} oand {Tn} of o and 7 respectively satisfy the moment equa- 


tion and [5.13] 
Remark 3.4 By deriving [3.9] with respect to x, and setting y'(r) — Z(x), one obtains 
the following second-order differential equation 
M [Z] (x) = la(x) Z" (x) + (x) + را‎ (2))Z' (x) = (Ansa - (0) Zr). (3.16) 
This leads us to denote and respectively by the weight and moment equations 
of M [Z] (x). 


Proposition 3.2 Let L|.| be the differential expression defined by [3.9] and @(.,.) is the 
bilinear form given in [.1] If L(p) = Ap and L(q) = vq for each p,q € P and À, v ER, 
with À # v, then رمام‎ q) = 0. 

For any solutions o and T of the Sobolev weight equations [3.11] and [3.12] respectively. 


Proof. This following from Theorem B.1] one has 


(A^—v)ó(p,q) = (An, q) — رمام‎ vq) 


= é(L(p), q) — é(p,l(q)) = 0. 
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1.€. 


(p,q) = 0. 


3.3 Classification of the Sequence Orthogonal Poly- 


nomial of Sobolev 


In this section, we consider the bilinear form $ defined by relation Knowing that © 
is a linear combination of two functional moments, we will try to classify the sequences 
of polynomials that are solution of the differential equation For this purpose, it is 
assumed initially that one of the two FM is quasi-definite and in a second time that the 


two are not quasi-definite. 


3.3.1 o quasi-definite 


Theorem 3.2 [49] Suppose that there exists an orthogonal SPOS (P,(x)],., with 
respect to @ defined by[3.1] solution of equation 3.9, If o is quasi-definite, then Ux) محم[‎ 


is necessarily classical with respect to cand T = kla(x)o for some real constant k. 
In addition, we have T = 0 if k = 0, or else T is also quasi-definite. 


Remark 3.5 [49] If o is quasi-definite and T = có, for all c € R* and for all a € R, 
then any OPS with respect to 0 given by B.1]is called S POS, but it does not satisfy the 
differential equation of Bochner's form [3.9] 

Moreover, this system satisfies another differential equation of order 2 with the coeffi- 


cient dependent on parameter n. 


Theorem 3.3 The classical orthogonal polynomials are also characterized by the fact 
that they are the only OPS satisfying the following condition: 
There exist 2 functional moments Tı not identically zero and To such that 


(Tı, PP.) + (To, PmPa) = 0, for m z n and n = 0,1. (3.17) 


3.3.2 T quasi-definite but c is not quasi-definite 


Theorem 3.4 [78 /9] Suppose that there exists an orthogonal SPOS (P,(x)),.g with 
respect to @ defined by [3.1] solution of equation [3.4 If T is quasi-definite then 

i) {Pi(2)}nsi i$ a COPS with respect to T verifies the differential equation [3.16] 

it) (P,(x)h,.y is a WOPS with respect to c. 

iti) lo(x)o = kr, for some k € R. Furthermore, either L(x)o = 0, or la(x)o. is 


quasi-definite. 
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3.3.3 c and ع‎ are not quasi-definite 


Theorem 3.5 Let {P,(x)},>0 be a sequence of orthogonal polynomials with re- 
spect to @ defined by [9.1] Ifo and 7 are not quasi-defined, then the polynomial sequence 
{Pn(t)}n>o does not satisfy the Bochner equation. 


3.4 The d—Orthogonality of Sobolev 


We give a generalization of the orthogonality of Sobolev, starting with the study of the 
regularity of the bilinear forms and then giving the conditions on the polynomial sequences 
d—orthogonal in the sense of Sobolev. We also show that the classical sequences are not 


classical in the sense of the generalization of the Sobolev orthogonality. 


Proposition 3.3 Let {P,(x)} 


Let V = (vo, v1,.….va 1) such as 


n>0 d— OPS with respect to regular U = (uo, ui, ….ug_1)T. 


Up = Ur + 0,6, , O<r<d—-1 (3.18) 


V regular if and only if 
1+ Lin+k-1(6, c) E 0 (3.19) 


where L\(x,c) is the polynomial defined as 


n—1 d—1 k 
Laine = Y, وه‎ 4 n RE PB) Let. (3.20) 
j=0 r=0 (Ur, Fate) r=0 (ur, Piel.) 


if o, = a then L5, (m, c) = aLan+k(x,c) where Las, (v, c) defined in 
Proof. Suppose that V is regular and let 


71-1 


Qm = Pm + > Omer Ami € R (3.21) 
i=0 
If we let m = dn +k with 0 > k <d-— 1, then we get 


(ur, Dane) E (ur, PantrPn) = ss for pm (3.22) 
í Qdn+r,dj+r (ur, Paj+rP;j) for 0 < J < n—1 
and by we also have 
(u Qa P;) = (Ur, O dein Pn UU a, Qan+r (6 PALO) for J =n (3 23) 
DR —aQan+r(c)Pr(c) fo O<j<n-1 
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Accordingly, we obtain 


Qan+n(Z) = Fantk(t) — Qanzk(€)Lanre-1(%; ©) (3-24) 
Set x = c in|3.24|to get 


Qaa (c) [1 + ايبوط‎ ©)] = Passo) (3.25) 
where necessarily 1+ L%,.,_1(c;c) # 0,otherwise we also have Pin+x(c) = 0. It now fol- 


lows by induction that c should be a common zero for mor than d consecutive polynomials 
{P,} wich is impossible ( corollary [2.3] ) = 


Let us consider d— bilinear forms رون ,... ,$1 روف‎ (d > 1). 
Consider a symmetric d— dimensional vector bilinear forms ¢ = (Qo, 1, ..., 04 1)7 
defined by 
Pr (p, q) = (or, p q) + apd), (3.26) 
where o = (09,01, ...,04-1) and T = (To, بك‎ ..., Ty_1) are d—dimensional functionals. 
The d—orthogonality in the sense of Sobolev is a generalization of the orthogonality in 
the sense of Sobolev. On the basis of the d—orthogonality definition, we give the following 


definition. 


Definition 3.6 We say that the sequence {R,(x)},-, is a sequence of d—orthogonal 
polynomials of Sobolev (d — SOPS in short) if it is orthogonal to ¢ which is defined by 


[3.26] i.e. 
Pal Rm, Rn) =0, n>md+a+1, m >0, 


(das Rm mate) z 0, m >0. 
We give necessary and sufficient conditions for the vector form © defined by to be 


regular. To this end, we study the following cases. 


(3.27) 


case 1: © = (20,01, ..., 04-1) is regular. In this case, let {P,},-, be the d—COPS with 
respect to c. 
We notice that if L,(r,y) is an orthogonal polynomial for some n > 0, then its 


derivative is given by 
NS dr*s 
LOU (x,y) = ——— Lr (x,y), > 0. 
ny) dx" dys (x,y), n2 


Theorem 3.6 The symmetric d— dimensional vector bilinear form ف‎ = (Do, $1, ..., $4.1) 
defined by 


$«(p,q) := (or, pq) + ى6)ه‎ pd) p,q €P and (<r > 0-1 (3.28) 
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is quasi-definite if and only if 
1 + aL?” (e, c) #0, m>1. (3.29) 
If ó is quasi-definite then 


0,1 / 
لكاي الم 


RACE) = Fats 
1+ aL (e, c) 


jor n>0. (3.30) 
Proof. For the first part , we assume that c is quasi definite. Then we can write 
m-—1 
Hj, = Pa + > gas (3.31) 
i=0 


If we assume that m = dn + k with 0 < k > d — 1, then we get 


T3 Page Pha f | = 
(or, Ponte) = K nd = b 5 )3.32( 
Qdn+r,dj+r (Or, Page) for 0 < J < n—1 
From , one has 
r(Ran+r, Pa) OR زعا‎ Pil for j=n 
(Or, RäntrP;) = ei n | anı (0) (e) ١ (3.33) 
0. (OPC) fo 0<j<n-1 
Accordingly, we obtain 
Ranse(®) = Passa) — aa, (o) Loss i5). (3.34) 


Differentiating and setting x = c, then we get 


Bins (6) [1+ LEP هايم‎ | = Pisa (9) 


where necessarily Tib cote c) # 0, otherwise we also have P7, (c) = 0. It now follows 


by induction that c should be a common zero for more than d consecutive polynomials 
UP as Which is impossible. Hence, follows immediately. m 

case 2: T = (To, Ti, ..., 74-1) is regular. In this case, let {Q,},, be the 0-025 with 
respect to to T. 1 

Now consider the following symmetric d—dimensional vector bilinear forms © = 


(Do, %1, ..., 94-1) defined by 


br (p,q) := (adc, p q) + (pq). O<r<d-1,ceER (3.35) 


Now we let {Qn(Z) محم[‎ be d-OPS with respect to 7 satisfying 


0 if n> dm +r +1 
kn(Tr) if n= dm +r, m > 0 


(Tr, Qn dm) E | 
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and {Ra(t)}h>0 be d- OPS with respect to Q satisfying 


0 if n>dm+r4+l1 


(Tr, RnR) = ; 0 
ka(dr) if n=dm+r (m 20 


The next proposition will prove useful later on. 


Proposition 3.4 The vector Q is quasi-definite if and only if a > 0 and 7 is quasi definite. 
Then, 

Rale) =0 for ner, 

kn+1(@) = (dm + r)mkam+1(T), 


1 
Qn = R! 


nti n+1° 


Proof. Assume that @ is quasi definite 


0 if n>dm+r4+1 m>0 


T Rm, Rn = 
Prk ) 8 if n= dm +r ,m >0. 
On the other hand 


(aô(x = c), Ii + (Ti Rao Ra) 


E | AR) Rs (c) + (ts RR) 


In particular, for m = 0 we have 


0 if n>r+1 


oa) = auto) = | OF 


l.e., 


R,(c)=0 for n>r+1 


Thus for m and n > 1 we have 


0 if > d 1 > 0 
Gin a) = Fl) - À au he eee 


kn(or) if n=dm+r m>0 ` 
Then, 7 is quasi definite and { R/,(x)} is d—OPS relative to 7. In this case we get 
Hop ue n>1 


and 


0 if n>dm+r+i 
(dm + r)Mkam+r(Tr) if n=dm+r,m >0 


(Tr, Rm Rn) = (TMNQm—-1Qn—1) = 
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Now, we assume that a 4 0, 7 is quasi definite and {Hn} >o is MPS, such that H, = 
nQj, 4. This leads to 


aif m=n =0 

0 if m=0 andn>r+l 

0 if n > dm +r +1 

(dm + r)mkam+r-1(Tr) if n= dm + r,m > 0 


br( Hm, Hn) = @Hm(C)An(c)+(t, HaHa) = 


ie. (H,),.g is a d- OPS with respect to © and by unicity, we obtain H, = Rn. This 
completes the proof. 
In this cas, we shall prove that the 2-classical orthogonal polynomials with respect to 


T, satisfy homogeneous third order linear differential equation of Boukhemis type. 


on (7) Quis (x) + analt) R(E) + anal) (2) oa(m)Quas(m) = 0, — (3.36) 
with o4 ;(r) for i = 1,..., 4 are polynomials of degree less than or equal to i and m 
an alz) = Fais 
is not define a 2—classical OPS in the sense of Sobolev 


Theorem 3.7 Assume that a 40 and T is 2-classic moment functional vector. 
If degani(x) > à — 1 then the sequence {R,(x)},-0 satisfies the following non- 


homogeneous third order linear differential equation 


Gy aR?) + (Qn-1,3 — alam "p CHEP + 04-12 — of, REF 


/ (3) y ES 
(an-11 — Ong CR pr 05, 4,3) n — Un, 


(3.37) 


where 


Un — Os LR s-- (Q1, — 07, wR +(e! sa din n > o (3.38) 


If degos;(x) = i then the sequence {R,(x)},>0 satisfies the following homogeneous 


four order linear differential equation: 
aaa + CT + Msn + 0511 43 = 0. (3.39) 


Proof. Let {Qn} >o be a classical orthogonal polynomial sequence with respect to 7 


satisfying and Ri (x) = nQn-1, n >1. 
For n = 0, Ro satisfies 
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For n > 1 and dega,;(x) € i — 1 one has 


(os i 4 RO) + (Qn-1,3 — l2) RI + (05 14 ie — an) 


= (n + 4)(o 14 QU), + (Qn-1,3 — Maes + (an-1,2 + CAM = "(مجم0 (و يرنه‎ 

= (n+ 4)(æ, 1409, + Gy sO + 051201, T (05, 12 + ars = Qr, 13) Qn42) 
3 n 
= (n + 4)(—an-1,1Qn-1 + (05, 15 + ara — 05 4,3) Qn42) 
= ( 


H 


/ (3) / 
—Q@n-1,1 Ter da O5, 13) Haa 


(3.40) 
and 
Gant + (An—1,3 — akaa la + (o 14 + Qn-1,2 — aR a 
= J (-an-11 + O12 + au = da او‎ 
= (Oy eu T a; = an 13) Ens + Un 
i.e., 


dnia RS + (Qn-1,3 — à ad) s + (0 14 F an-1,2 — d aR كي‎ 


/ (3) H u 
(an-11 — 0513 — Q4 F 05, 4,3) Rn+3 — Un 


for some constant vn. Since R,(c) = 0, n > 1 we have [3.38] 


If deg o4; (x) = i and according to we obtain [3.39] 
Thus the proof is complete. m 
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Characterizations of classical d-OPS 


4.1 Introduction 


Polynomials are special functions that we can differentiate and integrate many times. 
Any polynomial may be written as a linear combination of all other polynomials of lower 
degree, i.e., each polynomial of degree n must satisfy a linear recurrence of order n (n+ 1 


terms). If the recurrence is finite, i.e., there exists m > 0 such that 
n—2 
Pp (x) = (ant + bn) ER + 2 an kPk, 
k=m 


we say that the sequence of polynomials {P,,},, verifies an orthogonality of such type well 
defined. 

All notions of generalization of the orthogonality introduced by the authors, i.e., mul- 
tivariate orthogonality, multiple orthogonality, bi-orthogonality, d-orthogonality, ..., are 
obviously for searching a solution of some problems like Padé approximation, rationality 
or irrationality solution of higher order ordinary or partial differential equations, ... 

Classical orthogonal polynomials play an important role among orthogonal polynomi- 
als. The concept of d-orthogonality, d being a positive integer, appears as a particular 
case of the general multiple orthogonality. In particular, they satisfy orthogonality con- 
ditions with respect to several positive measures. Some recent works were focused on the 
analysis of properties of d-orthogonal polynomials generalizing the classical orthogonal 
polynomials and some characterization theorems were derived. 

Many examples of multiple orthogonal polynomials as well as of d-orthogonal polyno- 
mials were derived. They satisfy a linear differential equation of order d + 1, i.e., when 
d — 1 we meet the standard orthogonality. In this case, we know that there exist only 


four families of continuous and four families of discrete classical orthogonal polynomials 


43 
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solution of the Bochner’s equation, i.e., a linear equation of second order of the form 


az (x) y" + Bi (x) y + yy = 0, (4.1) 


where o» and f, are polynomials of degree two and one respectively and where ^ is 
constant. 

The solution of the equation are known as classical orthogonal polynomials. Start- 
ing from this point of view, many characterizations of the classical orthogonal polynomials 


are given. 


In this chapter, we wish to generalize most of these well-known characterizations for 
standard orthogonality in the case of d—orthogonality. Most of the results are given in 
the case d, but according to the complexity of the calculations we give some results in the 


case d — 2. 


Notice that the application of the classical polynomials is very large and in various 
domain. This is our wishes for this study. We hope that our results enlighten such a 
way for the applications in a new domain that not possible in the case d — 1, at least at 
the moment, we know that the classical 2 -orthogonal polynomials satisfy a third order 


differential equation which is impossible by the standard orthogonality. 


4.2 Preliminaries and notations 


4.2.1 m-Symmetric sequences 


Problems related to symmetrization of sequences of orthogonal polynomials on the real 


line play an important role [24]. In general, we have 
Definition 4.1 [16] Let m be a nonnegative integer. 
(a) A polynomials sequence {Pn }n>o is called m-symmetric if 


P,, (wet) = wy’ Pn (x), n 20, where wp = exp {2kri/ (m--1)), 1 <k € m. 


(b) A form T = (To, ..., Da 1)? is called m-symmetric if 
(D,, z 0*0) — 0, Qux m, vZp, n2 0, 
for each 0 € v € d — 1. In this case, 


(Ty, P (z)) = wy (uy, P (uwxx)), where wy = exp (2kri/ (m + 1)). 
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For the particular case: m = 1, we meet the well known notion of symmetric PS [24]. 
The components of the sequence {P,},,59 are (m +1) sequences {P/},59, 0 € p < m, 
such that 


Pm+1)n+u = 2" Pf peo) , 0 > ,روه > سر‎ n>0. (4.2) 
Definition 4.2 A d-OPS {P,},50 is called classical according to the Hahn property 
if Qn = (n+1) ! P}; is also a d-OPS. 


Since the derivative operator takes each d-symmetric PS into a d-symmetric PS, then 
a d-OPS {P,},50 is called d-symmetric classical, if Q, is d-symmetric d-OPS. 
A characteristic property of m-symmetric PS is given by the following 


Proposition 4.1 (16) For each d-OPS (Pj, with respect to U, the following state- 
ments are equivalent. 


(i) The form U is d-symmetric, 
(ii) The sequence {P,},50 is d-symmetric, 
(iti) The sequence {Pn},>o satisfying the recurrence 
Parasi (£) = ZPnra (£) — P1 Pn (a), n>0, (4.3) 
Pe wu 0 > 5 > 4. 


Note here that if {P,},>9 a m—symmetric d—OPS, all the components (P7],.o, 
0 X u € m, are d—orthogonal and if moreover {Pn},>, is classical, then all the (d+ 1) 
components are classical. 


4.2.2 Characterization of m— symmetric 


Let d be a positive integer and m de a nonnegative integer satisfying m < d. Next, 
we give a necessary condition on m and d to have an m— symmetric d—OPS and two 
characterizations of m— symmetric d—OPS. We denote by X* the multiplication operator 
by X^ in P 


Theorem 4.1 [16/Let {P,},59 be ad— OPS . Then the following properties are equivalent: 
1. The PS {Pn},>o is m—symmetric 
2. d +1 is a multiple of m + 1, say d +1 = p(m +1), and the PS {Pn},>o satisfies a 


(d + 1)—order recurrence relation of type 


p 
XP, = fad + S ng aa 


j=1 
with Yn p Z 0 and the convention P-n = 0 for all n > 1 
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Theorem 4.2 [16] Let {Pn}n>o be a d—OPS .Then its components {PF} ~o k = 0,..., m, 


n>0 " 
are d— orthogonal 
Classical d-OPS 
Theorem 4.3 [16] If {Pn},>o is a d— symmetric classical d— OPS, then its the compo- 


nents 23 AF k = 0,...,d, are classical d— orthogonal 


Theorem 4.4 [I6] Let {P,},>9 be an m—symmetric classical d— OPS then its first com- 
ponent {Pf },>o is a classical d- OPS 


4.2.3  d-quasi-orthogonality 


In the rest of this section, we recall the notion of d-quasi-orthogonality and some it 


characterizations. 


Definition 4.3 À sequence محم[‎ is said d-quasi-orthogonal of order s with respect 
to T = (To, sali) و‎ if for every 0 € a € d — 1, there exist s, > 0 and oa > Sa integers 
such that 

Ta, P4P,) 20, n>(m+s)d+a+il, m >0, 

(Ta, Pr, Pegs ata) À 0, mz 0, 
for every 0 < a > ل‎ — 1. We put s = maxs,, 0<a<d-1l. 


(4.4) 


Definition 4.4 [54|A sequence {Pa}n>0 is said strictly d— quasi orthogonal of order s 
with respect to I' = (To,…, Pedy if it satisfies 

(Ty, PnP 20, n>(m+s)d+v+1, m20, 

(Ta, PP) * 0, mz 0, 


for every 0 > v < d — 1 with s = Maxo<v<ad-1 Sa 


(4.5) 


Definition 4.5 Let v € P' and d, s > 1 integers. The sequence {P,}, is strictly 
s/d-orthogonal with respect to v if it fulfills: 


(v, PaPa) =0, nimd-s, m> 0, 
(4.6) 
(v, Pq Pmdars—1) z 0, m > 0, 

With this definition we can express: {Pn},s9 is d-OPS with respect to IT = 
(To, ql oa) if and only if is strictly (a +1) /d-orthogonal with respect to Ia, 0 < 
a €d- 1. 

A d -OPS محم[‎ with respect to V ,can be d—quasi orthogonal of order s with respect 


to — of dimension d . We have the following result. 
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Proposition 4.2 Let {Pa}nso 4 d -OPS with respect V . Suppose that exist T of 
dimension d such that ids is quasi orthogonal of ordre s with respect to T. Then 


necessarily d > d 


We consider only the case d — d , where we need the following characterization . 


Theorem 4.5 For any sequence {P,},50 a d -OPS with respect V , the following 


statements are equivalent 


(a) There exist T of dimension d such that UP aso is quasi -orthogonal of order s with 


respect to T 


(b) There exist T of dimension d and Sa, Ca > Sa integers, 0 > a < d — 1 such that 


(Ta, Pr) 20, n> sad+a+1, 


(4.7) 
(Ta, Poe Peas) A 0, 
for each 0 <a > 0- 1 
(c) There exist — of dimension d and d? polynomials 6%,0 < u, a > d—1 
0-1 
Li Diva, 0 > به‎ > 0- 1 (4.8) 
a=0 
where 
deg $ = Sa, 0 > بن‎ > 0-1 ond if d >2, 
deg كل‎ < Sq, 0<a<u-l, if l<a<d-1l (4.9) 


degó? > بن‎ - 1, utl<u<d-1 if 0<a<d-2, 


(d) There exist T of dimension d such that is strictly quasi orthogonal of order s with 
T 
respect to I = (To, ..., Fa_1) 


(e) There exist —of dimension d and sa > 0,integers, 0 < a > d — 1 such that 


Ta, Pr) 20, n2 sad + ديه‎ 1, 
To; P. dta) Æ 0, 


(4.10) 


for each0<a<d-—1 
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Lemma 4.1 A sequence {Pa }zosatisfy [4-4 can not be d—quasi-orthogonal of order 
s — 1 with respect to T = (To, al ea 


Definition 4.6 A sequence {P,,},,59 is said to be weakly d-orthogonal of index (k, 1) 
with respect to I' = (To, ..., Du if there exist integers, k,l > 1 such that 


Ta, P) 20, n> d(k,—1) ca 1, 

(4.11) 
(Ta, Patka-1)+a) * 0, 
(Ta, P) 20, n> d(la -1) +a +1, 


(Ta, X Pau, 3)4o) # 0, 


with k = max ką and 1!=maxl, for every 0 > a € d — 1. 


(4.12) 


Then, a strictly d-quasi-orthogonal sequence of order s with respect to I, is weakly 
d-orthogonal of index (s + 1, s + 2) with respect to U. 


Theorem 4.6 For each sequence {Pn} 


ing statements are equivalent. 


n>0 d-OPS with respect to —, then the follow- 


(i) There exists L € P' and an integer s > 1 such that 


GP, ح‎ 0, ns and (L,P;-1) 0, (4.13) 


(ii) There exist L € P', an integer s > 1 and d polynomials Q^, 0 > a < d—1 such that 
L= ncm tu, with the following properties: if s -l=qd+r, 0 <r <d—1, we 


have 
deg ¢" =q, 0<r<d-l, and if d> 2, 


deg $^ < q, 0<a<r—-1, if 1 > 7 > 0-1 (4.14) 
deg > ح بن‎ 1, r+1<a<d—1, if 0<r<d-2, 
Another characterization of d-quasi orthogonality in terms of sequences is the following 


Proposition 4.3 [70/For any two d—OPS,s {Pn},>o and {Qn}nso relative to T and Y 
respectively , the following are equivalent 


a T and V satisfy the relations [2.15] and [2.16] for some non zero d x d matrix polyno- 
mials © = (9) 


b there is nonnegative integer l such that 


dl 
P, (£) = Qui) + + Gens where ana > 0 (4.15) 


i=1 
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4.3 Classical d-orthogonal polynomials 


The characterizations of classical orthogonal polynomials in the algebraic aspect made by 
Maroni in [57], can be broaden here in the sens of d-orthogonality. We point out that 
the results given in the theorem below were the first characterization of d-classical OPS 
pioneered by Maroni and Douak (see, e.g., ). Here we first aim to simplify their proof 


and we shall give another characterization of d—orthogonality. 


Theorem 4.7 For any d-OPS {Pn},>g with respect to T = (To, ..., Ta)”, the next 


statements are equivalent: 
(i) The d-OPS {P,},=0 possesses the Hahn's property. 
(ii) There exist two d x d polynomials matrix © and للا‎ such that 
(ET) + VT = 0, (4.16) 


where ® = ($2) and V are in the following forms [30] 


0 1 0 
dex 0 
v= 4.17 
0 EE 0 d-1 
(rz) & cc Ga 


with w (x) is a polynomial of degree 1, [2 are constants and 


0(2) ° d$ (x) 


TIL m 


Poult) ° gaile) 


where $" are polynomials such that: 


degġ > 1, Usv<a+1, if 0<a<d—2, 
deg, =0, a+2<v<d-1, if 0<a<d—3, (4.19) 
degóÓ , <2 and degéót5 <1, 1> >ن‎ 0-1 


(iii) (P), sg is d-quasi-orthogonal of order 2 at most with respect to Y = ET and {Qn} 
is d-OPS with respect to V = ET. 
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Proof. we prove the implication (i) = (ii). Suppose that {P,},s9 is classical d-OPS 
and {Q,} d-OPS with respect to Y = (vo, v1, ..., v4 4)" where Q, = (n4- 1)! Pi ia BY 


theorem (2.11) we have 
Vdn+p = Ph Uo + Pn Ui کا‎ Po vai 


then, taking the derivative and using again the system (2.16) and using also (2.3) we 


obtain 
— (dn + u + 1) Pans, = (69) vo + ( ET tati lu 
— dé Po — De [Eria tdia] Pei 


Piet = Philo GLua rt os ile = a + alie + Galen 
(4.20) 


which implies 


(92) vo + (ph) vi +... + (68532) ون‎ = [dO — (dn + u + 1) 90,4] uot 
ie [C + Da + désigne —(dn+u+1) pi] Ut: 


note by 
ko(n,u) : = dé, Ÿ = (dnt p+ 1) Cas 
and 
Ring) : = ah? y + d£. — (dn + u عد‎ 1( eins 
with 


<n+1 if u=d—1 
and 
<n-1 if w+2<a<d—land 0<p<d-3 
deg[k.(n,u) = 4 <n if 1<a<pu+1 and0<y<d—2 
=n ifu=d-1land l<a<d-1l1 


If we take n = 1, u = 0 and (oi) = 1, we have 


d—1 


vo = ko(1,0)uo + V ^ ka(1, O)ua. (4.21) 


a=1 
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And for n = 1, u = 1, we obtain 
0-1 
(¢ ون"(‎ + v = Ko(1, L)uo + ذأ‎ Ko(1, ua. (4.22) 
a=1 
to simplify calculations take (01) = 1 for i € j 
1.e. 
0-1 
vı = [ko(1, 1) — ko(1,0)] uo + >} [ka(1, 1) — ka (1, 0)] ua 
a=1 (4.23) 
put n=1,1< u <d- 1, we obtain 
0-1 
v, = [ko(1, p) — ko(1, u — 1)] uo + 2 [ka(1, 4) — ka(1,  — 1)] ua (4.24) 


Note by 


ka(1,0) := 6 a := 0,...,0-1 


lka(1, 4) — ka(l, u — 1([ := ¢ a:=0,..,d-1and u= 1,...d— 1 


with 
deg of <1,if 0 > عدن > به‎ 1 and O<p<d-2 
deg of =0 fu+2<a<d—1land 0<p<d-3 
degóz , > 115 1<a<d—1and degol <2 
hence we have obtain 
peior. 
No, by taking n = 0, u = d — 1 in (4.20) and using also (2.3) we obtain 


dé = coli alsace Ca 103-4 
—d [vio + AD + AgEa + ... Ag iai], deg =1 


i.e 
; 1 
Vi = TE [-dV4 To = (co + dÀi)Ti d ccs (ca-2 + d)a-1)l 4-1] 
= — [Yro + El 1 + — +11 41] 
whence (4.16) and (4.17). 
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(ii) 2 (iii) Since Y = 9T, then 


0-1 
(Va, Pn) = 2, Tu PhPn) 20, n2 dt pl, 


حير 
(M Para) — (Ta; OS P) * 0,‏ 


for 0 <a > 0- 2 and 


(v4, Pa) = TOR + Xu (lu GP) = 0, nz 2d+p+1, 
(va—1, Poa) = (To, $9_1Podtn) À 0. 


Then there exists 0 € t € 2 such that {P,,} is d-quasi-orthogonal of order 2 — t with 
respect to Y = D—. 

(ii)— (i) Trivial since by hypothesis both {P,},59 and {Qh},>9 are d-OPS. m 

Unlike in the case d = 1, if {P,,} is d-classical (d > 2), we can't conclude from Pearson 
equation that {Qn} is also d-classical. Indeed, since {Qn} is d-OPS with respect to 
y = ŞT, then 


(DU) = PV + OY = PV + OVP := Vy جه‎ OV = vo. 


Then, the following problem is open: the derivative of classical d-OPS is again classical 
d-OPS? 
We try to give an answer to this question in the present section. 


Fairly, it seems as the first look that the answer is negative! In fact, we have an 
example p.86, Tableau II]. In spite of that the family (E) satisfies a second order 
linear differential equation, meanwhile this family isn’t classical at all. 

Our results show that the solution of Boukhemis and Zerouki differential equation [19] 
are the only 2 -classical OPS for which any derivative of that sequence is again 2—OPS. 


For this end, we set the following definition 


Definition 4.7 A sequence of d-orthogonal polynomials {P,,} is called very classical d- 
OPS if {P,,} as well as it's derivative of any order are classical d-OPS. 


So now, the question is: determine all the very classical d—OPS. 
On the other hand, note here also that we can extract a simple characterization from 


proposition [2.6] as follows (see [57]) 


Proposition 4.4 Let {Pa}„>o monic sequence verifying such that 
(a) Xn+d-10 = 0, 0 > به‎ > n= 1 and ومس(‎ SU, u> 0, 


(b) {Qn}nso is d-OPS with respect to V. 
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Then, the sequence Cab is d-classical OPS. 


According to the last assertion of Proposition since {Ph},50 is d-classical and 
Qn = (n+1)7' B,,, is also d-OPS, and from theorem [4.7] we know that {P,},50 is d- 
quasi-orthogonal of order two at most, then we have the following characterization of 
d-classical OPS 


Corollary 4.1 Let {Pn},>o be d-OPS with respect to T = (Vo, ..., 74-1). Then ies 


is d-classical if and only if there exists an integer p < 24 such that 


Bts Y wish Wig (4.25) 


v=n—p 
Note that when d = 1, this characterization is given in [53]. And for d = 2, this 
result proved first by Maroni in [58] using another approach. 
The next characterization is presented for the semi-classical case in [71], we give here 
a direct proof for the classical character and it can be also used to prove a similar result 
in the semi-classical case[71] . 


Theorem 4.8 [71] Let {Pa}n>o be d-OPS with respect to T = (Viens pal. Then 
{Pr}nso ts d-classical if and only if {Qn} nso is weakly d-orthogonal of index (1, 2). 
Proof. Suppose {Pn},>o d-classical, then {Q:},50 is d-OPS with respect to V = $—. 


Afterwards, we have 
(Un Qu =0, n> dnd he 1 


(Us; x C sna) = 0. 
In particular, respectively for m = 0 and m = 1, we have 


(ug, Qn) =0, n>u+1, 
(vu, Qu) # 0, 
and from the (d + 2) order recurrence of {Qn}, we have 
(Un, TQn) =0, n>d+u+l, 


(vu TQd+u) = 5 gi (Up, Qu) À 0. 

Comparing with the Definition [1.6] we conclude that (k,l) = (1,2). 

Conversely, suppose that {Q,} is weakly d-orthogonal of index (1,2) with respect 
to Y and we show that {P,} is d-classical with respect to D', i.e., we will prove that 
I satisfy Pearson equation (4.16). Let W = (wo, ..., Wa-1) such that W = V’, i.e. 
(Wa, T) = (uj, T), Vr € P. Then 


(Wa, Pa) = —n (Va, Qn_1) =0, nz a2, 
(Wa, Poti) = — (a + 1) (va, Qa) À 0. 
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Then, there exists a < ta < به‎ + 1 such that 


(Wa, Pn) =0, n>ta +1, 
(Wa, Pia) * 0. 


Hence, by theorem [4.6} there exist d polynomials wk, a € u > d — 1, such that 


d—1 m 
Wa = ) „=o Val a 


with, if ta = qad +ra, 0 € ro < d — 1, we have 


deg Vi? = qa, (<ra <d—1, and if d > 2, 
deg V4 < da, 0cutzre,-—1, if 1 <ra > 0-1 
deg YË > qa =1, ratl<p<d-1, if 0 €r, > 0-2 


Since a < ta < به‎ + 1, we distinguish two cases: 


(1) If0 > به‎ € d —2, whence a € t, > نه‎ +1 > d. Consequently, qa < 0. 


(2) Ifa 2d—1,thend — 1 € مخ‎ € d, whence qq. < 1. 


If qa-1 = 1, then rg_; = 0 and £4, = d, and if qu, = 0, then rg; > d — 1. Which 
prove that all polynomials are constant except V9 , = ax + b. 


Now, it remain to show that the matrix V is of the form (4.17). Indeed, from (2.3), we 


have 


vl, 1 0 => 0 p 
v! 0 2 0. : 1 
y = 0 = = . 1 0 
: . 0 : 
i^a 0 acm d Ta 


and from (2.15) and (2.16), we obtain 


d-1 
la = 910 (x) uo + 253 1o (2) Ty 


where $$, (x) = cx + d and where Qf, (x) = € for 1 < u € d — 1. Hence 


0 1 0 0 
vo 0 0 2 lo 
vi ١ 1 
. _ a T | 
" | 0 zm 0 d-—1 Tua 
Y (x) &r € Gat 
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Now, we let the form W = (wo,...,wa-1)’ such that W = (2V)', ie. (WaT) = 
((œva) ,7( , Vr € P. So, since {Qn} is weakly d-orthogonal of index (1,2) with respect 
to V, then 


ds dere (rue) 2) = -n (Va, TQn_1) 20, n>d+a+2, 
(Wa; Para+1) = — (d + a + 1) (Ua, x Qao) 7 0. 
Then, there exists نه‎ < t4 < d + a +1 such that 
(wa, Pr) =0, n>te +1, 
(wa, Fla) À 0, 


In addition, there exist d polynomials 0", a < u < d — 1, such that 


d—1 


Wa = < ر ر‎ 
put ta = qad + Ta, 0 € Ta < d — 1, we have 
deg 0? = qa, 0<r,; > 0-1 and if 0 < 2, 
deg 04 < a 0<u<ra—l, if 1 <ra < 0-1 
deg ð” «q4,—1, Tra+tl<u<d-1, if OSra <d—2. 
Since به‎ < t4 > d + بم‎ +1, we distinguish several cases: 


(1) 1150 > به‎ > 0- 2, whence a € t, <d+a+1<2d. Consequently, qa < 1. 


(2) If qa = 1, then ta = d+ ra > 0 +a +1 hence ra € ع به‎ 1. If qa = 0, then necessary 


Q X Ta. 
(3) Ifa =d -— 1, then d— 1 € t4-ı € 2d, whence qq_1 < 2. 


If qa-ı = 2, then rg, = 0, and if rg; < d — 1, then qa < 1. 
If رونو‎ = 0, then ra-ı = d — 1. This show that only 09 , is of degree two. 
Now, we have Y' = WI and (zV)' = OT. If we let V = $T, then 


O- = (VY = zV' + Y = UT + 9T. 


Hence 6 = © — z with Y = OF and (TY = VT. m 

Now for any sequence of OPS it's necessary to know their explicit expression of the 
measure for which the orthogonality is satisfied. It is however good enough to give a such 
tool for determining the orthogonalization measure. The quasi-monomiality principle is 
a very useful tool as well. Because it's not required any information about the classical 


character. Our approach deals only with the classical OPS and it will be possible to 
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generalized for semi classical case. In the following, we prove a new characterization of 
2—classical OPS, which will be without a doubt the pillar in this context. Note here that 
the next results are used before by Maroni and Douak to give an integral representation 


for some particular class studied starting from the corresponding generating functions. 


Theorem 4.9 For any 2-OPS (P, ], ss with respect to I = (Vo, T1)”, the next statements 


are equivalent: 


(i) The functional moment T satisfy the Pearson equation with (4.17), (i.e. V is 


classical], 


(ii) There exist fai}, and {8}; polynomials of degree n, and m, Ho, K2, Mo and 
N3, Lz of degree > 1, 2, 3 respectively, such that 


o (x) To + a3 (x) To + a5 (x) To = 0, 7. 
Bs (x) T1 + Bi (zx) Ti + B8 (x) T1 = 0, 
with 
Hı (x) r: = M (x) lo + N3 (x) Fg: 
and 


H, — w'ó —w(e —wv), Kə = whi — $n, La = w0 — OT, 
pa = Om, — 7), M» = m (d! — v) — 06, Na = Th — w0, 


Proof. Write the matrices in (4.16) when d — 2 in the form 


1 
Į = ? and 6 = n om) 
w(x) € (x) O(a) 
with € is constant and the polynomials v, w, m and 0 are of degree 1 at most and deg © < 2. 
In this case, the system (1.16) is equivalent to 


ul To 4- «o + 11م‎ + 214 = 0, 


(4.28) 
(9! — v) To + T6 + 41 + OF} = 0, 


with mı = 7 - 1 = C* and 0, 20 — .عأ = ع‎ 
First we establish all possible cases of the second equation in (4.26). So, eliminate T6 
from (4.28) we get 
Ho (x) To = Ko (x) li + L3 (x) T, (4.29) 
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where 


H3 — w'ó —w(d — Y), Kz = wh — pri, Ls = w0 — on, (4.30) 


with deg Hə < 2, deg Kə < 2 and deg L3 > 3. 
We distinguish three possible cases according to the degree of the polynomial Hp (x). 
1) deg Hp = 0, i.e. Hp (x) = C** :— C. In this case, replacing To and I, obtained from 


(4.29) in the second equation of (4.28) we find 


with 
D3(t) = 5 (9 — V) Ks + SOR - 6 
B, (x) = b (d — v) Ls + فط‎ (Kao + Li) + 0 (4.32) 
As (z) = GOL. 


Now it is necessary to search a common factor to reduce the degree of the polynomials 
As, B4 and Ds. In fact, we have 


Ds (x) = (x) 83 (x), Ba(x) = ¢ (x) 67 (z) and As (x) = (x) 8; (x) 


with 
84 = چ‎ [( — 20) m +w'6] (4.33) 
1 = Sl — 29) a+ Ko + 240 + wh — م‎ | 
Hence, we have (4.26) with deg 23 < 3, deg 87 > 2 and deg B1 < 1. 
2) deg Hy = 1, put H5 (x) = t. In this case, we have 
tTo = (HT 0)’ — HoT). 
Eliminating To from the system (4.28), we obtain 
then 
ipeum Ne OG He rie (4.35) 


Using again the second equation of the system (4.28) to get 
n (z) TY +n? (x) T; m (2) 1 = 0, 


where 
n? = (D - Y) [i (Ka — R1) Ho + 6m — w61] = (¢ — v) 83 
n? = (¢' — y) [I (Ko + L5 — Ro) Ha + on — w0] = (d! — v) 82 
7] 


1 —ll(g — yp) HL = (à —¥) 8? 
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with deg 63 < 4, deg 87 > 3 and deg (3 > 2. 


3) deg H = 2. In this case, we have by using (4.34) 
H}To = (Kli + LE) — Riri = ST 
Multiplying the first equation of the system (4.28) by H2H we obtain the second 
equation of (4.26) where 
Bi = ان‎ (d —w)(K,-— Ri) Ho + WHR) + Ti HH, 
pz = و اننا‎ = v) (Kə + L} = $3) H^» + wu H5» + T H2H} 
B3 — انا‎ Hala, 
with deg 82 > 5, deg 82 > 4 and deg 6} < 3. 
For the first equation of (4.26), eliminate w from from (4.28) we get 
where 
ln = 051; — 701, Ma =7 (¢ — v) — Ou", Na = دقع‎ w0, (4.37) 


with deg 1 < 1, deg Mə € 2 and deg N; > 3. 

We distinguish two possible cases according to the degree of the polynomial y (x). 

1) deg ui = 0, i.e. py (x) = C** :— C. As in the first part, we obtain after simplification 
the first equation of with 


aj = Tı 'ف)‎ —v) — w0 + M; 
o2 = om, — w0, + Mo + N} 
a3 = N3, 


where deg a? < 3, dega? < 2 and deg < 1. 
2) deg ıı = 1, ie. y (x) = k in the same way we have the first equation of (4.26) with 


od = TTE + LT + 7104 
a2 = WU + imu (M3 + N3 + ony = w01) +T (w0 = $71) 
o3 = tpa Na. 


where deg a? > 4, deg a? > 3 and deg < 2. 
Conversely, it is Hn to = that we can deduce the system (4.28) from (4.27). 


Indeed, on account of (4.29)-(4.30) and (4.36 4.36] )- (4.37), the system (4.27) may be written 
0 lwo + Ml: + rl] = wW (g = v) Do + 0414 + 0'1] ; 


(4.38) 
0 [mà + wTo + wro] =T [011 + (¢ = v) Do + orol š 
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After adding wéI and 015 respectively, to both sides of the first and second equality 


of (4.38), we obtain 


OX =wY 
OX = بلاج‎ 
with 


X -—[|wTs- ore + ml, + 211] 
Y = [(Ø' — p)To + T$ + €T. + #71]. 


Hence eliminating X and Y, respectively, we deduce 
(w0 — 16) Y = (w0 — xp) X = 0, 


whence the result since w0 — نع‎ # 0. m 


4.4 Application 


Furthermore, in [12], the authors showed that there are exactly 2° sets of d-symmetric حل‎ 
classical polynomials (see also [19]). The method used to obtain the differential equation 
of classical 2-OPS in [19], can be used again for the d-symmetric d-classical polynomials 


to give the explicit form of solutions. Indeed, we have the system 


= Yn loai A neas 


0 | 0 ~0 | 0 
in which we distinguish d + 1 cases in each case we have L solutions, hence we have 


E) 


solutions in total. In other words, solutions are 


So) Fou, = 49,, (trivial solution 
n+1l n+1 


0 — (+1+00 ~0 z0 — L0 

S1) Fanta = ( as, ) "Ydn-cag ANd Yanta = Vanta fora # ao and 1 < a < d 
~0 = n+l Pag 0 z 0 — n+1+Pa1 0 

52) Van+ao m ( N+ pog jte, and Yantay m ( n-Fpo4 Ydn+aı and 


Vata = Vina fOr a ينه روه عو‎ and 1 < a < d 


Sa) rg = ra for 1 <a < d. 


4.4. Application 


Riordan Arrays and d-Orthogonality 


5.1 Introduction 


5.2 Riordan arrays 


The Riordan group introduced by Shapiro et al [73] is a set of infinite lower-triangular 
integer matrices where each matrix is defined by a pairs of formal power series g(z) = 
sso gnz” and f(z) = SJ, faz” with go £ 0 and fi # 0. An infinite lower triangular 


is called a Riordan array if its i column generating function is 


matrix D = [dn,kln.k>0 
و‎ (x) [f (x)]’ for à > 0 (the first column being indexed by 0). With little loss of generality 
we also assume doo = go = 1. The matrix corresponding to the pairs g, f is denoted by 


(g, f). One example of a Riordan matrix is the Pascal matrix 


1 0 0 0 

1 2 1 1 0 0 
P= = 5.1 
(5) 1 2 1 0 


for which we have 


The group law is then given by 


(g,f).(h,l)=(g(ho f),lo f). 


The identity of this law is I = (1,x) and the inverse of (g, f) is (g, f) |! = (1/ (go f), f) 
where f is the reversion or compositional inverse of f. The reversion of f is the power 
series f such that (f o f) (x) =a. We shall sometimes write this as f = Revf. 
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A Riordan array of the form (g(x),x), where g(x) is the generating function of the 
sequence ay, is called the sequence array of the sequence ap. Its general term is a, ;. 
Such arrays are also called Appell arrays as they form the elements of the so-called Appell 
subgroup. 

For a Riordan matrix D = [dx], 459, Rogers [67] has found that every element dn+1441 
can be expressed as a linear combination of the elements in the preceding row starting 
from the preceding column. Merlini et al have found that every element in the 0 
column can be expressed as a linear combination of all the elements of the preceding row 
also see [75]. That is, there exist unique sequences A = (ao,a1,...) and Z = (20,21, ...) 
with مه‎ Æ 0, zo zz 0 such that 


(1) dn+1,k+1 = > don bis (k,n = 0, 1, و(‎ 


j=0 
oo 


(2) dn+1,0 = 3 glass (n = 0, 1, sl 


1-0 


The coefficients ao, a1,... and zo, 21,... appearing in (1) and (2) are called the A- 
sequence and Z-sequence of the Riordan matrix D = (g(z), f (z)), respectively. Letting 


A (z) and Z (z) be the generating functions of the corresponding sequences, we have 


We therefore deduce that 


and 
Z(z) =1/f (3) (1 - V/ (go f) ()): 
A consequence of this is the following result which was originally established by Luzón: 


Lemma 5.1 Let D — (g,f) be a Riordan array whose A-sequence, respectively Z- 
sequence have generating functions A(x) and Z(x). Then 


® 7 ةن دور 
DA AJ‏ 
If (ao, a4,...) is a column vector with generating function A (x), then multiplying‏ 


D = (g,f) on the right by this column vector yields a column vector with generating 


function g(x) A(f(x)). 
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5.3 Riordan arrays and d-orthogonal polynomials 


Let {Patnso be a sequence of monic polynomials in P for which there exist complex 


sequences {Bn}nsor (Xn), 0 € v € n such that 


Po(z)=1, P(x)=x- bo, 
n—1 (5.2) 
Patt (£) = (x — Bn4i) Pn (x) — 2 Xn (x), n 2 0. 


The dual sequence {Un}, >o: un € P' of {Pa}n>0 is defined by (un, Pm) := 0554, n, M > 0. 


In this case, we have 


Bn = (Uns E Pr (x)) >on > 0, 


7 (5.3) 
Xn» = (Uv, TPay1(x)), OS ou <n. 


Let us consider d forms uo, ..., رون‎ (d > 1). Let us now recall the following characteri- 


zation which we need in the sequel. 


Theorem 5.1 Let {Ph},>0 be a monic sequence of polynomials, then the following 


statements are equivalent. 
(a) The sequence {Pn},>o is d-OPS with respect to U = (uo, "e 


(b) The sequence {Pr},>9 satisfies a recurrence relation of order d + 1 (d > 1) : 


2-1 d—1—v 
Pai (x) = (£ — Bmta) Pm+a (£) — a كر ل ور‎ m+d-1-v (x), m < 0, (5.4) 
with the initial data 
P(x)=1, P(x) =z -— bo, 


5.5 
Py (£) = (£ — Bm-1) Pm-1 (£) — De, yee), 2E ms, Pu 


and the regularity conditions 49, , #0, m > 0. 


It is well known that we can express (5.4) as xP = طول‎ where P = (Py (x), P, (x) ,..,)* 
and Jq = (ai) 5-0 is a (d + 2)-banded lower Hessenberg semi-infinite matrix, i.e., 


for j >i+1 and i> j +d, (5.6) 


Qii+1 = 1, 1 > 0. 
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More precisely, the matrix ول‎ may be written 


Bo 1 Ü ^: 0 
TO 8 0 ee 0 
a 1° هم‎ 1 0 0 
Bal 4 HM Fe. ty ® (5.7) 
4 e ss 4^ quo À 0 
0 D + Vari Bari 1 0 
"r a 


and is called the monic Jacobi matrix of the monic d-OPS {Pn}, >o: 
Now, when {P,},,59 is not monic, then the linear recurrence (5.4) become 


Pn (x) and (ot LE 6n-1) Ph (x) = > yi n—l-v (x) Qn > 0. 


Write 
Prlz) = Does asdf ess (5.8) 


then 
ano = CnQn—1,0; 
"P (5.9) 
E 
üni fU d uisi = nait 1$ E 


On the other hand, for a family of monic d-OPS given by (5.4), and by using (5.8), we 


can have 


n+1 n n—1 n—d 
3 TONER = (x — Ba) 2 An px" = yet 3 s Lx =. Yay 3 And kt" 
k=0 k=0 k=0 k=0 
from which we deduce 
0-1 
An+10 = —Ênan,0 — D EE اوا‎ (5.10) 
k=0 
and 
d-1 
An+1,k = An,k-1 — na nd D WS Gus (5.11) 
k=0 
We note that if 6, and 4275, 1 < i < d are constant, equal to 8 and 74", 1 < i < d, re- 
spectively, then the sequence (1, —8, —4-1, —,4-?, ..., -7°,0,0,...) forms an A-sequence 


of the coefficient array. The question immediately arises as to the conditions under which 
a Riordan array (g, f) can be the coefficient array of a family of d-OPS. A partial answer 
is given by the following proposition. 
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Proposition 5.1 Every Riordan array of the form 


1 a 
5.12 
7 + »» brk 1 + xm x) ( ) 
is the coefficient array of a family of monic d-OPS. 
Proof. The array 0 (1 + YN ur à a (1 YN ur T) [40] has a C-sequence 
C(z) = > 5,59 nz” given by 


x x 
1+ 52 at 1—20(x) 


and thus 
d+1 


= — , Opx! 
k=1 


This means that the Riordan array 0 (1 + SCT ur Ph uu (1 + YS gag Ji is deter- 
mined by the fact that 


ik = ii + J tres for n,k = 0,1,2,... 
i>0 


where an,—1 = 0. In the case of (1/ (1 + Y er ae x/ (1 LOT our ‘)) we have 


d 


dn+1,k = An,k-1 — 1 6 1د‎ fin ede (5.13) 
v=0 


Working backwards, this now ensures that 


Pub) (x = 01) P, -5 bosai P ne v 


where P (£) = Sopp an, k£”. m 

We note that in this case the (d+2)-term recurrence coefficients B, and 4275, 1 > i < d, 
are constants. 

As an example, we have the following result 


Proposition 5.2 The Riordan array is the coefficient array of the modified 
d-orthogonal Chebyshev polynomials of the second kind given by 


PAV) = Un (£ — 01), n= 0,1,2,... 
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Proof. We have [31] (2.4)] 


1 oo 
= 27 t”. 
1 — xt + p» CRE AR 2. 
Thus 
1 OO 
= Un (x = 01) t” 
1 — (a - 01) + Ega Ong at 2. 
Now 


1 1 t 1 
1- (x - 004+ >= kytt 7 E + La TELE 1+ 3m iu) beum 


Thus 


0 3 = Un (x = 01) t". 
E + Yo et^ 1+ Eko n) Po 2 


a 
Also we have the following motivating result 


Proposition 5.3 Every Riordan array of the form 


d+1 
1— Ya Arf 3 )5.14( 
d+1 k? d+ k 0 
1+ مذ‎ One® 14-15, One 


is the coefficient array of a family of monic d-OPS. 


Proof. We have 
d+1 


B=(1 -X Ara", x) -A 
k=1 


where oe 
E m 1— nut Anc x 
B = (hs) = (Se aaa) 
A = (ang) = (37 F1 pk? 1+ THE 
and (1— 2 Axx^, x) is the array with elements 
1 0 0 0 0 0 
—AÀ 1 0 0 
— 2 —AÀ 1 0 
—À2 —AÀ 1 
—Agi : LINE e | oa 
0 —Àqd4i ۰٠۰ =۸ —A1 1 0 
0 


(5.15) 
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Now, by using the recurrence (5.13), we deduce 


d 
bn+1,k = dn+1,k — Y c جوم‎ ET 
d 
= An,k-1 — 2.28 Ov+1ün—v,k 


d 
—Ai 5 = pm NL 


d 
—... — نبول‎ |o. = 2 6 و0 1ن‎ ae " 
d 
= Dn,k-1 De Av+1bn—v,k- 


Hence, since bn, is expressed in terms of a,x, then as in proof of the previous propo- 
sition, B is a Riordan array of a family of d-OPS {Qn}, i.e. 


Qo(x) = 1, Q(x) = x — ĝi — À, Q»(x) = 2 — (20, + 3)z + A105 — ول‎ + 02 — Do, 


Proposition 5.4 The Riordan array is the coefficient array of the modified d- 
orthogonal Faber polynomials given by 


P,(z) = F(z — 01), n=0,1,2,... 


Proof. The d-orthogonal Faber polynomials {F, (x)},,59 are studied by Douak and Ma- 
roni [32], and in another way by Ben Romdhane [14]. This family of polynomials is defined 
with the help of the generating function as follows 


1 — yw Adr oo 
E =) Bf. 
l= at + Sopr kt 2 
Thus zu 
1— Ya Ad > 
aa E 
1= (g —61)t + aa Orat"? n=0 
i.e., 


l3 AQ t 1 = 
eh التشد‎ : = DJ E (2-0) t. 
1+ 53.904407 1-4 75 ont) 1-a مك‎ 
Also note that, we can express the d-Faber polynomials in terms of co-recursive d-OPS. 
In fact, Ben Cheikh and Ben Romdhane have shown that the c-Appell d-OPS are 


expressed in terms of d-orthogonal Chebyshev polynomials of second kind as follows 


d+1 


F, (z) = Un (z) + 9 Us i (x). 
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In this case, since Uf? (x) = U, (x), it follows that the Riordan array (5.14) is the 
coefficient array of the co-recursive d-orthogonal Chebyshev polynomials of second kind 
given by 


d+1 


P, (x) = Un (Œ — 01) + 3 Us (x — 61). 


k=1 


5.4 Stieltjes matrix 


Let L = (l, x) be Riordan and L be the matrix obtained from L by deleting the first 


n,k>0 
row, ie. lnk = .ريبما‎ For example, if J is the identity, then 
0 1 0 
" 00 1 
Iz 
0 D 0 1 


Observe that L = IL. There exists a unique matrix Sr = (Snk)n k>0 such that LS; = L. 
That is, 


lnk = » c - for n > 1. (5.16) 
i>0 

We call this matrix the Stieltjes matrix of L. For example, the Stieltjes matrix corre- 
sponding to the Pascal matrix (5.1) is 


1 
0 0 ub 1 


= ررم 


We note also that the Sz is unique, i.e., Sr = Sy € L= K. 


Theorem 5.2 If L = (g, f) is Riordan and Sr is (d + 2)-banded lower matrix as in 


(5.7), then 


(a) SL is in the following form 


B® 1 0 0 0 0 
4^ B i 0 0 
yo ae 8 À 0 
= qe ارم‎ 8 (5.17) 
n c 78 d 
0 2 ges „3-1 8 1 
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(b) f = نه‎ (1+ Bf +76 1f2 + ... Pf) and 9 = 
(1 - 2B -Ef — ..— rf). 


Proof. Write S; in the form (5.7) with C; (x) the generating function of the itè column 
of L, i > 0. We have: C; = gf’. By looking at the first column of LS; and L, we obtain 


and 
Bo =ho=n 


liobo + : du = læ = go 
0-1 ded _ 7 
loo Bo + layi +9 = I39 = 93 


laobo + lat + laa- + 7? = laa = Japi- 
Multiplying the first equation by x, the second one by z?,.. and the last by x%+1 and 


summing we get 
zbog + ayî gf + ... + gf =g- 1. 


Whence و‎ (x) = (1 - zbo - هبه‎ if. zf) 
Repeat the procedure for i > 1, and in each case we try to find a system that contain 


all coefficient 6,771, ..., 49. One can see that we have 

GP = ral EIU exa. (5.18) 
that is 

Ge (Cii + BiC; + Cite hs + YP, Cira) . (5.19) 


Hence after simplification we find the expression of f. 
For instance, rewrite the equality (5.18) in terms of j and subtract, we get 


0-1 


(Bi = B;) f E ` e = a) for - 0 


v=0 


which is equivalent to 


Bi = Bj = B, quic suem 0 > >ن‎ 0-1 


0 


لا 
This leads to the important corollary which is an analogue of [81].‏ 


Corollary 5.1 If L = (g(x), f(x)) is a Riordan array and P = Sy is (d + 2)-banded 
lower given by (5.17), then L^! is the coefficient array of the family of d-OPS. 
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Proof. By Favard's theorem, it suffices to show that L^! defines a family of polynomials 
{P,(x)} that obey a (d+2)-term recurrence. Now L is lower-triangular and so L~! is the 


coefficient array of a family of d-OPS P,,(x), where 


1 P(x) 
z P(x) 
Lr r? = P(x) 
g” P3(x) 


We have 
رل د ررم‎ Sl shal Sb athe bh el si 
Thus 
Spek ele...) = LI. Tew) =L was 
We therefore obtain 
Bo 1 0 0 
a" B 1 0 P(x) zPy(x) 
de qe ug 0 P(x) P 
LO gy». yd- P(x) |= | xP(x) |, 
79 n, yl 1 P(x) xP3(x) 
0 4o . . 


from which we infer that 


Pitt) م ع‎ - bo, 
P(t) = (z — b) P, (x) — yf! 


Pia (£) = (x — 8) Pa (x) - YF Pee ay P (نه)‎ — Y8. 


and 
d—1 
Poste) = (x m 8)Pa(z) m $ 45 tas) 4 > 1. 
v=0 
لا‎ 


Combining these result with previous proposition, we have 


Theorem 5.3 A Riordan array L = (g(x), f(x)) is the inverse of the coefficient 
array of a family of d-OPS if and only if its Stieltjes matrix Sr is (d + 2)-banded lower 


Matriz. 
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We give the following example 


Proposition 5.5 The Stieltjes matrix of the inverse Riordan array 


1 x 
E Toast 1-45 i) 


left-multiplied by the k^ binomial array 


1 x = 1 x 
1—kz'1—kr] Ml-zrz1-zm 


61+ k 1 0 0 
0 A+k 1 0 
05 0 +k 1 


is given by 


Proof. We have 


1 7 1 z ) 
1+ Y One®’ 1+ DT Opa 1+kx 1 + تدا‎ 


1 x 
i (; + (01 +k) x + Sf Pera 1+ (0, + k) + iy اج‎ | 


And more generally, 
[1-511 Ag 2 ( 1 g ) 
San rot) te 1+ 


7 ( 1 — ART x ) 
N14 (0 E) z Y hna 1+ (0, +b) + Y a J 
Therefore, the inverse of the last matrix has the Stieltjes array 
6+ بل -يم‎ 1 0 0 
059—259 +k 1 0 
03 — ول‎ b +k 1 
03 02 


(5.20) 


0411 mm Ad+1 
0 04.1 


0 
= 
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5.5 Exponential Riordan arrays 


The exponential Riordan group is a set of infinite lower-triangular integer matrices, where 


each matrix is defined by a pair of generating functions g(r) = go + git + gos +--+ and 
f(x) = fix + fom +... where go £0 and fı Æ 0. In what follows, we shall assume 
go = fi = 1. 


The associated matrix is the matrix whose i-th column has exponential generating function 
g(x) f (x)!/1! (the first column being indexed by 0). The matrix corresponding to the pair 
f,g is denoted by |g, f]. The group law is given by 


l9, f] ١ [h. 1] = [g(h o f), to f]. 


The identity for this law is I = [1, z] and the inverse of |g, f] is [g, f] ^ = [1/(g o f), f] 
where f is the compositional inverse of f. We use the notation eR to denote this group. 
If M is the matrix |g, f], and u = (u,)1>0 is an integer sequence with exponential 
generating function U (x), then the sequence Mu has exponential generating function 
g(x)U(f(x)). Thus the row sums of the array |g, f] have exponential generating function 
given by g(x)e/(? since the sequence 1,1, 1,... has exponential generating function e”. 
We will use the following [25] [26], important result concerning matrices that are pro- 


duction matrices for exponential Riordan arrays. 


Proposition 5.6 Let A = (ax), >o = [g(x), f (x)] be an exponential Riordan array and 
let 


cy) = Co + رون‎ ey He r(y)m rot ry roy t+... (5.21) 


be two formal power series that 


r(f(x) = f(x) (5.22) 
_ g(x) 
c(f(z)) = 2" (5.23) 
Then 
(i) An+10 = 2 ciana (5.24) 
(ii) Antik = rod +o Dil cia + Erie) ang (5.25) 


i>k 
or, assuming Ck = 0 fork < 0 and rg = 0 fork < 0, 
1 f 
An+1,k = m 3 i!(Ci—k + kr;_x11)Qns. (5.26) 
i>k-1 
Conversely, starting from the sequences defined by , the infinite array (ak), 450 
defined by is an exponential Riordan array. 
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A consequence of this proposition is that the production matrix P = (Dis محر‎ for an 
exponential Riordan array obtained as in the proposition satisfies 


i! | 
Pij = ge +jriju)  (ea=0). 
Furthermore, the bivariate exponential function 


E xz 
Qp(t, z) m » Pn "i 


of the matrix P is given by 
p(t, z) = e" (e(z) + tr(z)). 


Note in particular that we have 


r(z) = f'(f(z)). (5.27) 
me (F(z) 
_ IMM 
c(x) = GF (5.28) 


Now, an important property of the matrix L which has led to the proof of various 


results in this context is the following 


- d 
Proposition 5.7 [//] We have L = TE 
x 


d 29 n pe n 


n=0 


Proof. Since 


then equating the first columns of the matrices L and LS, we obtain 


bo = ما‎ 
lobo +41 = lao 
lobo + la | + 7? = موا‎ (5.29) 


laobo + lave + ...laa-171 +9} = lau, 
and for n > d+ 1, we see that one can has 
loo + جيرا‎ + Inari + Hlina = .ممما‎ (5.30) 


Multiplying the second equation of (5.29) by x, the third one by z?,... and the last by z^ 
and (5.30) by x”, and summing we get 


Bo (1 E lox + loz? à sa) | 1 (x klar? 4 us] | EX (z^ + laprast +...) = g' (x), 
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Bog (x) -- o 9 (x) f (2) +... +729 (x) OF =o, 


Now, by equating the second columns of the same matrices we get 


L el 

ho + 81 = l3 

loo + labi +41 = la 
d d2 (5.31) 

lao + 13161 + وجووا‎ + 6 = la 

laxi + lata 81 + lass + o + laa +7! = lapon, 
and for n > d +2, we see that one can has 

lno + ln, 101 + bane +... + lna + lna) = سما‎ (5.32) 


Multiplying the second equation of (5.31) by x, the third one by z?,... and the last by z^ 
and (5.32) by x”, and summing we get 


(1 | hogz | Lot? | 8) | By (x | lar? | ل‎ 


Tg Gr oe) و‎ sas.) 


= [f + liz? + 4) ; 


g (x) + Big (a) f (x) +49 1g (x) [f (2)? + ... +989 (2) [f (1 = g (2) f (x), 
whence 
1+ Bif (x) +8 [f (a) + ...+ 9 [f GOL = f(a). 
لا‎ 


Theorem 5.4 If L = (g(x), f(x)] is an exponential Riordan array and Sy, is (d + 2)- 


banded lower matrix, then necessarily 


a 1 0 0 
ud; ay 0 0 
وم‎ M a 1 0 0 
SL = | "WI uw s Su 5 (5.33) 
م‎ Hy wy Ga 1 0 
0 qui py Hé” am À 
0 
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d—i 
where {ax}x>0 is an arithmetic sequence with common difference a, and [e ; 
> i) k>1 


1 <i > 0, is an arithmetic sequence with common difference 7? and 


In(g) = "nc tug f(x) + ... + mo lS Go) dx, 9(0) = 1, (5.34) 


where (a), is the Pochhammer symbol and f is given by 
d 
f =1+af +> Y P", f(0-0, (5.35) 
i=1 
and vice-versa. 


^ 1 
Proof. We consider the lower triangular matrix L with iij (x) [f (z)]" for the exponential 


generating functions of the kt? column, k > 0. We note that L is a Riordan matrix with 


exponential generating functions. 


From (5.35), we see that 
d 
[z"] g' = bole] و‎ + 9 92 ^ [2"] af", 
i 


In41,0 = Bolno + YF a + ... + Yin a: 


Hence, for k > 1, we have 


1 f 1, 1 m 
(ur) = gol gy t 
d - d o. 
- روق)‎ + Soft) oft + qy (hte Df) of 


ATEN 
then 
5 1 d 5 1 0 5 (B + kf) 
| ( di 00 - 1 0 
epo] EEE) pu .ل‎ par EE PD محرو‎ 
— n 1 | k—1 n Ok k | 
[x^] (ki + lej sof 
0-1 
TH sg pedis ارو‎ 
(+1)! (kd)! 7 
where 
ax = o + kB, and uL! = (k + 1, (+k °) for 1<i<d, 
whence 


ln+1,k = 506 + Qkln,k + u n1 Feat Ulo ktd- (5.36) 
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Finally, we easy observe that 


nt ome 
Qk — Ok 1— Û := and k — kl adi for 1 <i < d. 
لساك‎ EED (G ? Eu 


= 
We give now an example for a d-OPS defined by a generating function which can be 


also generated by a Riordan array. 
Let us consider à family of Hermite polynomials generated by the following generating 


function 
d 


G (x,t) = exp fi —b)t— 3 e E pM TP, (x) t”. (5.37) 


i=l 


This generating function satisfies the following first order differential equation 


G, (x, t) -(: —b— Yu Jet (5.38) 


Now replacing G (x,t) and ©“ (x,t) by the right hand side of (5.37), we can easily obtain 


the linear recurrence satisfied by the sequence of polynomials {Pa}n>0 as follows 


Parayı (£) = (x — b) Pasa ( 3T n +d); Para (x), n20 (5.39) 


which means that {Pn}, 59 is 2-OPS where (n +d);=(n+d—-i+1), 
By translation, we can take a = 0. In this case, we have G7, (x,t) = tG (x,t), i.e., 


Pi, (x) = (n + 1) Pa (£), (5.40) 


that is, (P5), is Appell sequence, and then it is classical d-OPS. 
Now by choosing 
b=d and a =i+1 for 1<i<d, 


the corresponding Stieltjes matrix (5.7) may be written 


b 1 0 ee 0 
01 0 1 0 ee 0 
2a» 201 0 1 0 -< 0 
6a» 
d'aq vis da, b 1 0 
0 ag(2+d), --- d(d+l)as (d+1)a b 1 
0 . . l ; à 
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where the common difference is 0 in all case, i.e., 
HR ME L 
(k+d), (k+d-1), 


Therefore, from (5.34) and (5.35) we obtain 


f(x) =a and In(g) = Je + at + 2agx”... + dlagr*)dz, 


a; — a; = 0. 


d 
g(x) = exp L + Yaaa ! 


i=1 
Corollary 5.2 {4JIf L = [g(x), f(x)] is an exponential Riordan array and Sr, in the form 
(5.7), then L~' is the coefficient array of the family of monic d-OPS. 


Clearly, L^! is then the coefficient array of the family of polynomials {P,(x)}. Gath- 


ering these results, we have 


Theorem 5.5 An exponential Riordan array L = [g(x), f(x)] is the inverse of the 


coefficient array of a family of d-orthogonal polynomials if and only if its production matrix 
P = S; is in the form (2.33). 
Proposition 5.8 Let L = [g(x), f(x)] be an exponential Riordan array with (d + 2)- 
banded lower Stieltjes matrix Sr. Then 
n d—1 a n 

EOS) — un") بح‎ 
where U = (uo, "m is the linear functional that defines the corresponding family of 
d-OPS. 
Proof. Let L — (les eso: We have 


Applying U, we get for each 0 < a € d — 1, 
Ual”) = ua (Dio ln; Pi()) = Yo lius (Pi) 
= Vico lnidia = lna = n![z"]g(z) [f E)" . 
= 
Corollary 5.3 [44/Let L = [g(x), f(x)] be an exponential Riordan array with (d + 2)- 


banded lower Stieltjes matriz Sr. Then the moments un of the associated family of d-OPS 
are given by the terms of the first column of L. 
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5.6 Sheffer Riordan array 


Let و‎ (x) and f (x) defined as above. Then the polynomials P, (x), n = 0, 1, ... defined 
by the generating function 
g (z) e — ^ P. (z)t", 
n20 

are called Sheffer type polynomials with P) (x) = 1. 

Therefore, the set of all Sheffer type polynomials (P, (x) = [z^] g(x) e*/? Y with an 
operation of umbral composition, forms a group called the Sheffer group. 

In the following paragraph, we give an example of Sheffer type Riordan array. Let 
us now consider a family of d-OPS of Charlier type defined by the following generating 
function 


n20 n! " 
i=1 


G (x,t) = (1 + rt) 9/? exp pe] =). Lg (x) t". (5.41) 


This generating function is of Sheffer A-type zero. Note also that when 7 — 0, we get 
the previous example. 
The family of polynomials generated by (5.41) satisfy the following recurrence relation 


Pau (x) = (x +b ح‎ Tn) Pa (x) — ` ai (n); Pai (£), n>0 (5.42) 
i=1 
which means that {P,},59 is d-OPS where a; = (n + 1) 0:41 +6; and (n); = (n — i + 1),. 


In this case, the corresponding Stieltjes matrix take the following form 


b 1 0 0 
01 b—T 1 0 0 
و20‎ 201 b— 2r 1 0 
6a3 6a» 301 b— 37 
d'aq dai b— dr 
0 ي(0 + 2)وه‎ d(d--1)a9 (d+1)a 
0 : 


d—i 
where the common difference of {ax = b},.9 is 7 and it is 0 for { E } i 
= k>0 


EA, 
Therefore, from (5.34) and (5.35) we obtain 

f(x) "كمس دمع‎ and 

hug) e [o + a4 (b + ze”) + 2a; (b + xe™) +... + d!aa (b + xe dx, 
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i.e., 


d 
lng (x) = 5 ixt exp {iTr}. 
i=1 


5.6. Sheffer Riordan array 


Conclusions and future directions 


In this thesis we are concerned with Riordan arrays and d-orthogonal polynomials. After 
giving some preliminaries results in the two first Chapters we focused in Chapter 3, 
on the construction of the inner products which generate the sequence of d-orthogonal 
polynomials in the sense of Sobolev. Our perspective in the future is to study d-classical 
of Sobolev and d-semi classical orthogonal polynomials of Sobolev. 

In chapter 4, we characterize the d— classical orthogonal polynomials. . We present an 
algebraic theory of classical d—orthogonal polynomials and we want to fill in some gaps. 
We broaden and close some inclusions that exist and are known perhaps as consequences. 
In the future we aim to extend the combinatorial theory of orthogonal polynomials created 
by Vienno to d-orthogonal polynomials. Indeed, a second method of calculating moments 
has been discovered by using a non-trivial combinatorial model from a generalization of 
Pollaczek polynomials that can be generalized to all Appell polynomials. 

In Chapter 5, we pioneer the study of d-orthogonal polynomials and Riordan arrays. 
More precisely, we have extended the results of Riordan arrays and orthogonal polynomials 
in such a way to give a significant generalization from the orthogonal case to the d- 
orthogonal case. We show that the Riordan array, in the constant coefficient case, is a 
good vehicle for such studies. A principal discovery is that the d-orthogonal polynomials 
are characterized by the Rordan arrays. Our perspectives are then to extend this idea to 
algebraic structures in d-Hankel matrices and d- Hankel-plus- Toeplitz matrices relating 
to d-classical orthogonal polynomials. We aim also to extend our research to d—continued 


fraction and d-semi classical orthogonal polynomials. 
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